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Abstract 


The  various  types  of  unified  field  theory  of  gravitation  and 
electromagnetism  are  reviewed,  and  that  which  is  presently  known  about 
Bonnor-Moffat-Boal  Theory  is  outlined  in  some  detail.  It  is  shown  that 
there  is  no  simple  method  to  generalize  the  concept  of  duality  of  the 
electromagnetic  field  to  nonsymmetric  theories  in  a  meaningful  way. 

Several  families  of  possible  solutions  important  for  analysing  the 
viability  of  Bonnor-Moffat-Boal  theory  are  described. 

An  attempt  is  made  to  physically  interpret  the  geometric  structure 
of  Bonnor-Moffat-Boal  theory  (and  by  generalization,  other  theories)  to 
the  extent  that  general  relativity's  geometric  structure  has  found  inter¬ 
pretation.  Many  difficulties  threatening  the  viability  of  the  theory  are 
found  and  left  unresolved. 

A  weak  perturbation  expansion  formalism  is  developed  and  with  it  an 
approximate  magnetic  monopole  solution  is  found  in  the  theory.  The 
geometric  optics  approximation  is  developed  and  applied  to  the  theory, 
to  study  the  propagation  of  light  and  gravitational  radiation.  In 
relatively  simple  situations  radiation  is  found  to  travel  along  null 
geodesics.  In  more  complicated  circumstances,  difficulties  are  encountered, 
and  no  conclusions  can  immediately  be  drawn  regarding  propagation  of 
radiation  in  these  cases . 
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Chapter  I 

Review  of  Unified  Field  Theories 


1 


2 


Introduction 

This  thesis  is  mainly  concerned  with  gaining  a  deeper  understanding 
of  Bonnor-Moffat-Boal  unified  field  theory  (BMB  theory,  hereafter)  of 
gravitation  and  electromagnetism.  In  the  latter  parts  of  this  chapter 
we  will  review  other  unified  field  theories  to  gain  an  overall  view  of 
the  subject  and  to  be  able  to  compare  and  contrast  BMB  theory  with  its 
competitors,  which  we  will  do  in  subsequent  chapters.  In  chapter  2  we 
will  go  on  to  review  BMB  theory  as  it  is  presently  known,  and  consider 
some  special  circumstances  in  which  new  solutions  of  the  equations  might 
(be  naively  expected  to)  exist.  In  chapter  3  we  will  consider  just  how 
much  of  general  relativity's  deep  geometric  structure  can  be  found  in 
unified  field  theories  with  .'no asymmetric  metrics,  BMB  theory  in  par-  ^ 
ticular.  In  chapter  4  we  will  develop  the  formalism  of  weak  perturbation 
expansions  about  exact  solutions,  primarily  as  a  step  toward  invest 
tigating  "geometric  optics"  in  chapter  5>  but  with  some  obvious  applica-** 
tions  of  its  own.  In  chapter  5  we  will  investigate  light  and  gravita¬ 
tional  radiation  within  the  framework  of  the  geometric  optics  approxima¬ 
tion,  mainly  for  the  purpose  of  deciding  whether  such  radiation  travels 
on  null  geodesics  or  not.  In  chapter  6  we  draw  such  conclusions  as  we 
can  from  the  work  done . 

Weyl's  Theory 

Einstein's  general  relativity  (which  we  will  also  call  Einstein-Maxwell 
Theory  and  abbreviate  EMT,  when  we  are  considering  gravity  plus 
electromagnetism)  couches  all  quantities  relating  to  the  gravitational 
field  in  terms  of  geometric  objects  on  manifolds.  The  gravitational 
effects  of  other  fields  are  included  by  identifying  their  stress- 


energy-momentum  tensor 


T  ,  with  the  Einstein  tensor  G  ,  : 
ab  ao 


(i.i) 


, 

■ 

•  . 

u 
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is  a  properly  geometric  gravitational  quant it lyf  but  T  has  no 
a  priori  geometric  content  whatsoever.  Thus  the  other  fields  do  not 
enter  in  a  geometric  way. 

The  only  other  field  known  at  the  time  general  relativity  was  first 

proposed  (it  is  still  the  best  understood  field)  was  the  electromagnetic 

field.  Several  researchers  including  Einstein  soon  began  to  look 

for  a  way  to  combine  gravitation  and  electromagnetism  in  a  "natural" 

geometric  way,  so  that  both  emerge-  from  geometry  on  manifolds.  Further 

motivation  for  this  was  provided  when  it  was  found  that  exact  solutions 

in  EMT  consistently  display  singularities .locations  where  the  curvature 

quantities  diverge.  It  was  thought  that  a  correct  theory  should  generate 

elementary  particles  as  regular  solutions,  and  that  modifying  EMT  to  a 

correct  "unified  field  theory"  might  get  rid  of  the  singularities.  So 

(1) 

far,  this  has  been  a  false  hope. 

The  first  such  unified  field  theory  was  Weyl's  (1918) ^  ^  . 

(For  review  of  Weyl's  theory  see  refs.  4,5).  In  it  he  proposed  what 
today  would  be  called  a  gauge  theory  of  length  on  a  manifold  with 
metric.  In  general  relativity  the  length  of  a  vector  is  preserved  in 
parallel  transport.  That  is,  the  inner  product  of  two  vectors: 

(1.2) 


Ta.'  T  a  .  b 

k  Ja  =  sabk  J 

of  which  length  is  a  special  case: 

-.2  .a.  _  -a  .b 

1  =  J  Ja  -  sab°  J  ’ 


(1.3) 


is  not  changed  when  we  move  the  vector  from  point  to  point  according  to: 

,  .a  r,a  .b,  c  (i  l) 

Si 

In  more  modern  notation,  for  transport  in  the  direction  of  k  ,  this  is 
expressed  as : 


t  b . 

k  J 


a;b 


0 


(1.5) 


. 
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(semicolon  denotes  covariant  differentiation  with  respect  to  the  connection 
).  This  is  so  because  we  arrange  that  the  metric  be  compatible  with 
the  connection: 


5ab;  c 


=  0 


Weyl  suggested 


(6) 


(1.6) 

letting  the  length  vary  with  displacement ,  and  as  a 


first  relaxation  of  invariant  length  we  let  the  variation  be  linear  in 


length  for  infinitesmal  displacement: 
dl  =  (I  dxa)l . 


a. 


(1.7) 


Note  that  the  new  vector  quantity  j  enters  (1.7)  in  the  same  manner 


that  the  connection  enters  (1.4). 


From: 
2 


a F=  21^  a  dxa  =  d(gabjaj1:’) 


we  can  quickly  show: 

.a  ,-a 


r^o  =rtc(1c)  -«  K +  gac  K-  ^ 

where  ra  \  is  the  standard  Levi-Givita  connection  found  in  EMT . 
bc(lc ) 


(1.8) 


(1.9) 


a 


Note  that  both  g  ^and  /  are  needed  to  determine Tbc 

Now,  consider  conformally  rescaling  the  metric  with  some  scalar 
field  f (xa) : 

sab  =  fSab  (f>0)  •  (1-10) 

This  modifies  lengths: 

A?  2 

lc  =  fl 


as  if  we  had  a  new  f  : 

O/ 

A  if 

<p  —  d  4-^  .a, 

'  o  'a 


A  Q 

rl  =  r 


a 

a 


a  =  /  ^  +2  (inf)  j 


a 


(l.n) 

(1.12) 

(1.13) 


be  '  DC 

(where  the  comma  refers  to  coordinate  differentiation) .  The  Weyl 
geometry  thus  admits  a  "gauge"  transformation,  and  we  call  /  a  the  gauge 
vector  field.  Quantities  which  do  not  change  under  such  transformations, 


' 

. 


5 


gauge  invariants,  will  be  considered  of  greater  physical  significance 
than  gauge  dependent  quantities . 

If  f  is  zero,  we  regain  length  preservation,  and  are  in  the  realm 

(?) 

of  general  relativity  again'  '  .  We  can  reduce  </>  to  zero  by  a  gauge 

3/ 

transformation  if  it  is  the  gradient  of  a  scalar,  which  is  equivalent  to: 
Fab^a,b-^,a  =  0' 

F  is  a  gauge  invariant  antisymmetric  tensor  of  the  second  rank.  By 

cL  D  ' 

construction  it  satisfies: 

F  ,  -  =  0.  (1.14) 

(ab ,  cj 

The  brackets  here  indicate  antisymmetrization  over  all  three  indices. 


kbk 


a;b 


We  can  still  identify  geodesics  using  parallel  transport: 

=  0  (1.16) 

but  these  geodesics  are  net  extremals  of  arc  length,  since  nonzero  arc 

lengths  are  now  gauge  dependent.  Null  vectors  and  geodesics,  however, 

are  still  invariant.  We  can  form  curvature  tensors  in  the  usual  way. 

As  an  example,  the  curvature  scalar  has  the  form: 

R  =  R  V  -  6  f  <j>  +  6^  (/-g1  ^a)„. 
emt  jzzp  ,a 

(R  is  the  curvature  scalar  as  would  be  found  in  EMT,  g  is  the  determinant 
emt 

of  the  metric) . 


3/ 

Let  us  label  an  entity  by  the  way  it  gauge  transforms.  If  T  *  _ 


transforms  as: 

A 


then  T 


ij,a .  .  . 

,3/  •  •  « 


b.  .  . 


(1.18) 

Is  said  to  have  a  weight  of  n  under  the  gauge  transformation. 


_  •  •  • 

b .  .  .  b .  .  . 


Note  that  raising  and  lowering  indices  changes  weight.  Some  examples  of 
these  that  will  concern  us  are  given  in  Table  1.  In  it  : 


Fab  =/=?Fab 


(1.19) 


The  underline  indicates  a  tensor  density . 


' 


’ 
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TABLE  1 


Object  Weight 


/=? 


..ab 


F' 


ab 


R 


bed 


R' 


ab 


2 

-2 

0 

0 

0 


R 


-1 
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A  conventional  way  to  find  field  equations  is  to  use  a  variational 


principle: 


I  4/W/=gd4x  =  0 


(1.20) 


(gab  and  wil1  be  varied  independently) ,  in  which  W should  be  a 
gauge  invariant  scalar  density.  On  purely  formal  grounds: 


+  wabggab)d4 


implies : 


wa  =  0 


x  =  0 


wab  =  0 


(1.21) 


(1.22) 


These  are  fourteen  equations,  but  they  do  not  comprise  an  independent  set. 
We  can  generate  identities  by  formal  gauge  and  coordinate  variations 
in  (1.20).  An  infinitesmal  gauge  transformation  implies: 

Wa  =  2Wa  .  (1.23) 

Infinitesmal  coordinate  transfomations  generate  a  further  four  identities. 
Thus  (1.22)  really  represents  nine  independent  field  equations. 

The  theory  begins  to  show  weakness  when  we  look  for  the  particular 
form  of  W.  It  is  reasonable  to  ask  that  W  be  constructed  from  the 
metric,  its  first  and  second  derivatives,  the  gauge  vector  and  its 
first  derivatives .  These  restrictions  give  us  only  four  rational 
functions: 


F  ,  Fab  ,  R  ,Rabcd  ,  R  ,Rab  ,  R2  .  (1.24) 

ab  abed  ab 

It  was  hoped  at  the  time  of  proposal  that  the  correct  unified  field 
theory  would  have  its  Lagrangian  uniquely  determined  with  no  unknown 
factors  (the  values  of  which  must  be  guessed)  and  with  the  gravitational 
and  electromagnetic  terms  appearing  together,  not  artificially  separated. 
Here  we  have  the  undesirable  situation  of  four  expressions  to  be 
combined  somehow,  including  one  (F  )  which  is  strictly  electromagnetic, 
Further,  there  is  no  way  to  construct  a  generalization  of  the  Lagrangian 


of  EMT: 


■ 


. 
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R  +  4  F  F' 
ab 


,ab 


(1.25) 

because  now  R  is  not  gauge  invariant .  The  closest  we  can  get  to  this 


is: 


w  =  -4  R2  +  A  F  ,  Fab 

4  ab 


(1.26) 

and  this  is  therefore  the  expression  considered  most  likely  to  succeed. 


An  immediate  problem  is:  this  will,  as  it  stands,  generate  gravitational 

field  equations  of  fourth  order,  instead  of  EET's  second  order.  Weyl 

(8) 

tried  to  get  around  this  by  a  particular  choice  of  gauge.  ' 

0  =  SJw/lgA  =  +  A^F^VgOlaV  (1.27) 

Since  R  is  of  weight  one  (  and  in  cases  for  which  R  is  of  single  sign) 


we  can  choose  a  gauge  so  that : 
R  =  -1  =  constant  . 


(1.28) 


This  is  a  measure  of  the  overall  curvature  of  the  spacetime,  somewhat 
in  the  sense  of  the  cosmological  constant .  We  can  then  bring  our  varied 


action  into  a  more  recognizable  form: 

_ab 


S([r  .  +  AF  ,Fab  -  3(?f  ^ X  =  0  . 
,u  emt  ab  4  a  J 


(1.29) 

The  normalization  used  implies  cosmic-sized  units.  To  go  to  conventional 


units : 


^a 

x 


€«1 


(1.30) 


with  simultaneous  rescaling: 
1 


^ab  £  2  ^ab 


(1.31) 

,  Fa'>  =  Fab  (1.32) 

R,mt  +  AFabFat  +  |2(1  -  ^/a)J/l?d4x  =  0  .  (1.33) 

This  looks  like  the  vacuum EMT  expression  with  correction  terms  on 


the  end  which  we  hope  are  small  because  of  the  smallness  of  the  cosmological 
term  e .  To  generate  Einstein's  cosmological  constants 


(1.34) 


p 
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we  set: 


X=  §f  . 

2 

By  varying  f  we  get  electromagnetic  equations: 

Si 


(1.35) 


t^ab 


=  s 


= 

A 


(1.36) 


and  charge  conservation  is  automatic : 

=  0  •  (1.37) 

However,  to  get  a  gauge  invariant  expression  for  (1.36)  we  use  the  gauge 
invariant : 


/v/a 

s 

-  .a 


4+1  (lnH)  a 


implying : 


't  =  zi.  034  +  jr  j  g aDy=g’ . 


ab 


>by 


Variation  of  the  metric  components  gives: 
R  ,  -  (R  +  X)  g  ,  =  AT  , 

ab  — - —  ab  ab 


‘ab 


=  (FodFCd  +  Sab  '  FacFbC  -  4Sb 


(1.38) 

(1-39) 

(l.ijfl) 

(1.41) 


(Q) 

Einstein  immediately  raised  an  objection,  based  upon  what  we  will 
call  "geometric  interpretation",  that  is,  based  upon  the  demand  that 
all  the  physically  interpretable  geometric  quantities  in  EMT  should 
have  analogies  in  unified  field  theory.  In  EMT,  the  ajec  length  of  a 
timelike  path  is  related  to  the  elapsed  proper  time,  and  the  length  of 
the  tangent  to  the  path  is  the  rate  at  which  the  proper  time  passes. 

If  we  apply  this  directly  to  Weyl's  unified  field  theory,  we  find 
that  the  rate  at  which  a  particular  clock  ticks  depends  upon  the  value 
of  {  at  the  clock's  position,  and  in  general  continually  changes.  We 

Si 


can  easily  show 


(10) 


then  that  two  identical  clocks  at  the  same  location 


with  different  past  histories  have  different  rates  of  ticking,  the 
difference  increasing  with  time .  This  past  dependent  behavior  is 
surprising  and  is  considered  unlikely.  An  interesting  reinterpretation 


. 
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was  provided  "by  Londorr  ' ,  who  showed  that  if  we  ask  for  geodesics  along 
which  tfiis  length  variation  does  not  occur  (thus  avoiding  the  above  stated 
problem),  we  start  to  get  quantum  mechanical  results,  eg.  Bohr’s 
quantization  in  a  spherically  symmetric  field.  However,  this  idea  was 
not  extended  far. 

While  the  generalization  of  EMT  that  Weyl  proposed  is  geometrically 
elegant,  the  Lagrangian  for  the  field  equations  is  not  uniquely  defined, 
and  the  resulting  field  equations  do  not  seem  to  be  a  simple  generalization 
of  those  found  in  EMT.  Possibly  because  of  the  complexity  of  the  field 
equations,  Weyl's  theory  was  soon  abandoned. 

Five  Dimensional  Theories 

The  five  dimensional  (including  "  projective")  unified  field  theories 
form  a  family  unrelated  to  Weyl’s  theory.  This  type  was  first  proposed 

by  Kaluza^^”'^  ,  in  1921,  and  later  extended  by  Einstein  and  collaborat- 

(1^-17)  i  (4) 

orsv  J  .  In  the  forthcoming  analysis,  we  will  follow  Bergmann  s 

review. 

These  theories  use  five  dimensional  Riemannian  geometry  in  which 
four  dimensions  are  distinguished  from  the  fifth.  We  will  begin  our 
review  by  considering  a  general  "four  dimensions  within  five"  formalism. 

In  a  five  dimensional  space  with  coordinates  y  ,  metric  g  ("<*”  in¬ 
dicates  five-space  indices  0, . . .4; "a"  indicates  four-space  indices  0, ... ,3) > 
we  introduce  four  parameters  x  ,  functions  of  y  ,  and  assume  their 


derivatives : 

Q,a  =  7.^) 

are  linearly  independent.  We  have  a  set  of  curves,  x  -  const,  which  fills 
-the  five-space,  and  which  forms  a  structure  on  the  space  invariant  under 
the  "parameter  transformations": 


r*  a  „a/  b« 

x 


=  fa(xU)  . 


(1.43) 


* 


' 


p 
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We  can  find  a  large  number  of  geometric  objects  in  this  formalism. 

There  are  ordinary  five  dimensional  tensors  (indicated  by  Greek  indices), 
and  p-tensors,  or  parameter  tensors  (indicated  by  Latin  indices),  trans¬ 
forming  as  four-dimensional  tensors  with  respect  to  parameter'transformations . 
Entities  with  some  Latin  and  some  Greek  indices  mix  the  transfomation 
properties  above .  We  identify  a  fundamental  vectorf ield  A*  (tangent  to 
the  curves)  by: 

(1.44) 

(1.45) 

and  define  a  "reciprocal  Q" : 

(1.46) 

Ak  =  0  .  (1.4?) 

With  these  we  can  decompose  any  vector  into  a  p-vector  and  a  scalar: 

(1.48) 


toe  A*  =  0 
A“A^  =  1 


a  c  a 
^=sb 

cx 


ka  =  <£  4* 


k  =  k 


o< 


and  conversely  associate  a  vector  with  any  p-vector: 


<x  „  <*  a 
m  =  Q  m  . 

cl 


(1.49) 


(1.50) 


(1.51) 


We  now  form  the  p-metric: 

Sab  Qa  SoP  ' 

The  mixed  Q-quantities  project  from  five  dimensions  into  four,  orthogonal 

CX 

to  A  .  The  projection  tensor  is: 

h  ^  -  Ql, 


=  Si-  V*- 


(1.52) 


Thus: 


k*  =  Q“ka  +  A“t. 
a 


(1-53) 

We  have  several  new  types  of  derivative.  First,  there  is  the 


"A-derivative  of  a  p-tensor": 


3/  •  •  • 


b . .  . , 


(1.54) 
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•  •  • 


If  the  A-derivative  of  k  ’ *  *  vanishes  it  is  "A -cylindrical"  (the 

U  •  •  • 

term  will  be  more  meaningful  later) .  Then: 

=  <Sr(Ari<r+  U-55) 

Fab=Qf^VV  (1‘56) 

5-  = A*V  -  V  (1'57) 

Fab,«A°<'=  ‘  (1'5'3) 

There  is  the  "p-derivative" ,  a  generalized  derivative  with  respect  to 

a 
x  : 

ex 


(1-59) 

If  k  is  A-cylindrical,  then  this  is  just  the  ordinary  four-space  coor  - 


k,  =  k  Q, 

J  a  ,cx  a 


dinate  derivative  ..  In  general  p-derivatives  do  not  commute: 

k.  ,  -  k.,  =  k  (Q'j  -  Qf:  ) 

jab  |  ba  ,  <r  a| b  b [a 

=  k  A^  F 
,<r  ba 

There  is  the  co variant derivative  of  a  p-tensor: 


.a  T  a  i-|U  T  c 

k  =  k  +  r  k 
;<r  ,  ^  1  co- 


(1.60) 


(1.61) 


Notice  the  mixed  indices.  In  specifying  P  it  can  be  shown  that  the 

a  <r~ 


choice : 


r^=  i 9  a  ^icy  ' 


satisfies  both  of: 


g  ,  =0 

ab;<* 

^o<  ^b ;  cr=  °  ’ 


(1.62) 

(1.63) 

(1.64) 


Meanwhile,  the  "covariant  p-derivative"  is  the  ordinary  covariant 


derivative  multiplied  by  Q  : 

cl 

T  cx  T  ,  p*  <r 

k  ;b  ■  k  |b  +  Hrb  k  ! 


(1.65) 


' 
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,o( 


-  r*  nf 


crb 


r;A 


.  a 
k 


=  ka,,  +  r~  k 


;b  |b 

K*  =  < 

_  i  ad, 

=  TS  (g 


a  -f  c 
cb 


(1.66) 

(1.67) 


ed|b  +  gdb|c  ”  Scb|d)  * 


(1.68) 


Because  there  are  several  derivatives,  there  are  several  curvatures: 

(1.69) 


T  C*  0<  c<  S 

k  V*  -ki^=R^k 


k  -  ka  =  Ra  ;'>a 

;/?y  dTj? 

R  =  Q"  (R  U/3  +  A  "  A^Ai;«') 

k°%=  -  k%  =  +  (Av  ■  V}  A"kVJ 

aS1 


R' 


L-a  _  Da  -*T  d  T  a 

k  ,  -  k  ,  =  R.  -  k  +  k  , 

;bc  ;  cb  deb  fd 


cb  * 


(1-70) 

(1.71) 

(1.72) 

(1.73) 


’^c^and  R  are  the  standard  four  and  five  dimensional  Riemann 
tensors,  expressible  in  the  usual  way  in  terms  of  the  respective 
connections: 


Rf 


) . 


deb 


.  ex'  . 
+  A  .A 


.n.p  —  A  A  c 
*\/3  \p 

For  Kaluza's  theory  we  will  use  the  double  contraction  of  this: 


(1.74) 


c*a  bc^d 

5  ds  R 


eba 

cr  „Kc*,ko- 


=  R  -  A  -  (A^  )2 

6  r;(T  K  \<r*J 


(1.75) 


+  B<rB--2k-(A^|CC+2fi^t 

R  is  the  five  dimensional  curvature  scalar  here. 

Most  of  the  theories  single  out  "special"  coordinate  systems: 
those  coordinate  systems  where  the  fifth  dimension  is  not  in  evidence, 
These  are  identified  by  the  demands  that: 


a  a 
y  =  x 


(1.76) 
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A*  =  (0,0, 0,0,1). 

Transformations  between  special  coordinates  are  of  the  form: 

T  =  ^V) 

y5  =  r1  +  f5(yb) 


coupled  with: 


r\jQ. 

X 


=  fa(xb) 


In  special  coordinates  we  have; 

Wa  ,  1) 

v =  h* +  -4 

^  -i 


g 


ab 

g 

be  / 
-g  P, 


ac  / 

-  -  8  ?c 

-  i+«°Vd 


(1.77) 

(1.78) 

(1.79) 

(1.80) 

(1.81) 

(1.82) 

(1.83) 


(1.84) 


/  transforms  according  to: 

K =  ^  -  f5.b} 

so  that  f^  induces  a  "gauge"  transformation  on  (  (but  not  in  the  same 

cl 

sense  as  Weyl's). 

In  limiting  our  attention  to  special  coordinates,  the  multi-faceted 
"four  dimensions  in  five"  structure  we  have  just  discussed  becomes 
considerably  restricted.  "A-differentiation"  becomes  differentiation  with 
respect  to  y^: 


a.. .  .  a  **  _  v3,  •  •  • 

^  "I  **  *1  r-* 

b...}of  D.  .  .  O 


For  p-differentiation: 

a  _  c  a 
Qoc  - 

q?  =  sr  ,«*=  i.2, 3.^ 

cl  d 


,«=  5 


(1.85) 

(1 .86) 
(1.8?) 
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k,  =  k  -  <L  k 
la  ,a  ra  ,5 

(1.88) 

Then: 

Fab  =i4,b  ‘  4,a  +  44,5  '4 

>a  ,5 

(1.89) 

B*  =  W  W(5,  0) 

(1.90) 

Art;«=  r5=<  =  ln  (  det  V)  } 

*5 

(1.91) 

and  in  fact: 

iet  =  det  gab 

(1.92) 

A°<;«=  *(ln  (det  8ab)  )  >5 
i  ab 

(1.93) 

=  *g  gab,5  ' 

Now,  Kaluza's  Theory  assumes  that  the  five-dimensional  metric  is 
A-cylindrical  -  a  strong  assumption.  It  implies  that  is  a  Killing 

cL 

vector,  vanishes,  and  that  the  A-curves  (the  original  x  =  constant 

curves)  are  geodesics.  F  ^  is  A-cylindrical  and  satisfies: 


F  =0 

[ah,c3 

and,  in  special  coordinates: 


(1.94) 


Fab  ^a ,  b  ^b ,  a 

To  get  field  equations  we  assume  a  variational  principle  with  the 
Lagrangian: 

X  =<1  ldet  6^1'  E 

in  analogy  with  general  relativity  (gravitation  only).  Now: 

cba  +  ^ 


(1.95) 


„  <?a  be  _d 

R  =o  dS  R 


ab 


(1.96) 


(1.9?) 


In  special  coordinates  we  can  replace  det  g  by  det  g^,  and  because 

the  Lagrangian  is  A-cylindrical,  the  integral  may  be  taken  over  a 

a 

five-volume  or  over  the  four-volume  domain  of  the  parameters  x  : 

(sbcRaota  +  y  ,yb )  =  °  (i-98) 


S/(, 


ab; 


where : 


'  ,  ' 

* 
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g  =  det  gat  . 

(1.99) 

The  variation  is  subject  to: 

(Ssab)  =  o 

(1.100) 

(g^a>,5  =  °  ' 

(1.101) 

This  variational  principle  and  the  resulting  field  equations  are  identical 
to  those  of  EMT. 

The  projective  theories  of  Veblen  and  Hoffmann^  ^  and  of  Pauli 

are  reinterpretations  of  Kaluza's  theory:  these  authors  really  only 

worked  in  different  types. of  coordinates.  The  A-curves  are  assumed  to 

represent  only  one  point  of  the  physical  space .  Then'  the  metric  must 

he  A-cylindrical,  and  the  field  equations  are  again  equivalent  to  EMT. 

Since  Kaluza's  theory  is  equivalent  to  EMT,  it  does  not  predict 

anything  new  (and  does  not  solve  the  problem  of  singularities,  for 

example).  We  can  of  course  expect  new  phenomena  in  modifying  it. 

(15) 

Einstein  and  Mayer  tried  introducing  a  five  dimensional  formalism 
in  four-space  without  actually  introducing  five-space.  Then  the 
"five-dimensional  vectors"  transform  according  to  arbitary  matrices 


M  : 

Y°'=  (1.102) 

If  these  were  true  coordinates  transformations  M6^,  would  be  restricted 
by  some  differential  identities.  This  construction  creates  a  large 
number  of  new  differential  covariants,  and  it  is  difficult  to  attach 
meaning  to  all  the  new  quantities. 

Another  interpretation  of  the  formalism  appears  in  theories  in 
which  the  fifth  dimension  is  closed.  These  were  proposed  by 
Einstein,  Bergmann,  and  Bargmann  starting  in  19381  •  To  constrict 

a  space  with  closed  fifth  dimension',  we  cut  a  thin  slice  out  of  a 


* 


. 
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five-dimensional  oo ntimiion  and  identify  the  two  open  (four-dimensional) 
faces,  resulting  in  a  "tube".  We  assume  that  closed  geodesics  that 
go  around  the  tube  intersect  themselves  smoothly. 

The  term  "thin"  applied  to  the  slice  means:  thin  enough  that  there  is 
exactly  one  closed  geodesic  in  the  fifth  dimension  through  each  point . 

The  arc  length,  S,  of  one  circuit  of  the  closed  geodesic,  the 
"circumference"  of  the  tube,  is  a  constant  of  the  space.  The  vectors 
tangent  to  these  closed  geodesics  form  the  unit  vector  field  A  : 

A*  A*  =  0  .  (1.103) 

;  o 

In  special  coordinates  this  becomes: 


-  =  0  .  (1.104) 

a  >  j 


3 

The  rest  of  the  metric  components  are  periodic  in  y  ,  and  indeed  any 
field  uniquely  defined  in  the  space  is  periodic  in  y^  with  period  S.  Then: 

=  0  (1.105) 


(1.106) 


l  ,  -  =  A  -A 

ab,5  a; b  b;a 


(1.107) 


If  we  wish  to  derive  field  equations  from  a  variational  principle, 


we  have  four  scalars  of  second  differential  order  to  choose  from  (  all 
others  differ  from  these  by  divergences)  to  construct  the  Iagrangian: 


~a  be  _  d 

$  d  S  Ralo 

(1.108) 

Fab  F  , 
ab 

(1.10.9) 

(A  +  A  )  (A  c  +  Ar  )  g  ^ 

p\°< 

(1.110) 

(A  .  +  A  )  (A  ^  +  Aj  )  g**  g^  . 

v  «\/s;  /?;<* 

(1.111) 

In  the  variation: 

o 

ii 

!>i 

ux 

Ti 

qjs 

(1.112) 

(this  is  stated  in  special  coordinates,  with  g  =  det  g^ 


)  we  must 
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preserve  the  geometry  of  the  closed  space.  In  special  coordinates, 
must  he  independent  of  y-5,  and  must  he  periodic  in  y-5.  These 
restrictions  clash  with  the  usual  requirement  that  the  variations  vanish 
on  the  boundary  of  the  domain  of  integration.  Instead,  we  can  ask  that 
the  boundary  extend  exactly  once  around  the  tube,  and  that  the  variations 
vanish  on  the  part  of  the  boundary  generated  by  the  A-curves.  This 


.  modified  procedure  unfortunately  results  in  integro -differential  equations: 

/(PabSgab  +  Ja£Va)  m?  a5y  =  o  (1.113) 

* 

which  implies: 


(1.114) 


s 

f  Ja/=?dy5  =  0  . 

5 

y  -0 


These  satisfy  five  integro-differential  identites: 
s 

rb\  r~~\  ,  5  _ 


(2  Pa  ;b  +  KJ  =  ° 


ba 


yJ=o 


(1.115) 


(1.116) 


/  (jb;b  +  ,5)  /^dy5  =  0  .  (l.ll?) 

y-^o 

—  3j 

P  anc^  ^  generally  contain  the  expressions  in  the  field  equations  of  EMT 
(among  other  things) . 

The  appearance  of  integro-differential  equations  can  be  objected  to 
on  the  grounds  that  physical  fields  are  usually  described  by  pure 
differential  equations.  With  a  closed  fifth  dimension,  such  equations 
probably  cannot  be  found  in  a  straightforward  way  from  a  variational 
principle . 

Einstein's  Theory  of  The  Nonsymmetric  Field 

The  most  popular  unified  field  theory  is  Einstein' s^ 20 ^  (1945) "Theory 
of  the  Nonsymmetric  Field".  In  it,  one  still  has  a  four  dimensional 
.  manifold,  but  the  metric  is  no  longer  symmetric: 


f 
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Sab  S(ab)  +  s 


[ab] 


g(ab)  =  i  (gab  +  gba}  * 
—  ~2  _  )  • 


(1.118) 

(1.119) 

(1.120) 


g[abJ  2  ^&ab  &ba 
The  antisymmetric  part  is  expected  to  be  related  to  the  electromagnetic 

field.  The  line  element  in  the  theory  is  expected  to  be: 


ds^  =  g  ,  dxadx^ 
Dab 

,  a,  b 
dx  dx 


(1.121) 


=  O' 

o 


(ab) 


In  the  original  proposal,  Einstein  raised  the  possibility  that  g 


(ab) 


is  real  while  g^-j  is  Imaginary,  so  that  g^  then  has  "hermitian 


symmetry”: 


Sab  =  gba 


(1.122) 

where  the  bar  indicates  complex  conjugation.-  '  The  interpretation • ^  „ 


generally  favoured,  however,  is  that  g  ^  is  strictly  real  and  this  is  the 
interpretation  we  will  use  unless  otherwise  stated.  Since  the  BMB  theory 
is  a  modification  of  Einstein's,  having  different  field  equations  in 
the  same  general  framework,  some  discussion  relevant  to  Einstein's 
theory  will  be  found  in  later  chapters.  In  chapter  3>  we  will  dis¬ 
cuss  some  of  the  interpretational  problems  of  the  nonsymmetric  theories 
in  general,  including  the  particular  problems  of  the  hermitian  symmetric 
alternative . 

We  assume  a  nonsymmetric  connection  : 
rf  =Fft  S  +  frf  -J  (1.123) 

1  be  (dc;  {be/ 

where  f^a  \  transforms  as  a  connection,  but  P as  a  three  index  tensor, 

the  "torsion”.  From  the  connection  we  find  the  curvature  tensor: 

a  _  pa  _  pa  +  na  p®  -  Pa  pf  (1.124) 

bed  bc,d  bd,c  'ed'bc  1  gc1  bd 

and  the  contracted  (on  upper  and  last  lower  index)  curvature  tensor, 


no  longer  symmetric: 


' 


' 
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R  ,  —  R/  ,  \  +  R  _ 
ab  (ab)  [ab] 


_  p ^  _  n ^  _j_  pC  p b  _  pG  pd 

ab,c  '  ac,b  '  ad  '  cd  '  ab  /  cb 
Before  we  can  identify  the  curvature  scalar,  we  must  specify  the  upper 

indexed  metric.  We  have  two  choices: 

(1.126) 


(1.125) 


ac 


s  ®bc=Sb 


S  S0b  =  Sb  ■  (1.12?) 

both  choices  involving  the  assumption  that  det(g  )  is  not  zero, 

3j  D 

reasonable  if  the  metric  is  to  raise  and  lower  indices.  The  former  alter¬ 
native  is  usually  preferred  r(also  -  it  generates  the  proper  EMT  limit 


in  BMB  theory)  and  will  be  used  in  the  rest  of  the  thesis.  Then: 


R  =  gab  R 


ab  • 


(1.128) 


For  field  equations,  Einstein  chose  his  Lagrangian  in  analogy  with  the 


S/b 


B  Tpg1  d^x  =  0  . 


pure  gravitational  field: 

(1.129) 

cl 

We  use  B  for  the  curvature  and  also  for  the  connection  in  order  to 
distinguish  them  fro'm  transformed  quantities  which  will  be  used  in  the 
field  equations.  We  do  not  assume  any  a  priori  relation  between  g  ^  and 

O  ^ 

L^c,  but  vary  them  independently.  The  interdependence  of  g^  and 
is  one  of  the  field  equations.  Performing  the  variation,  we  find: 

sab,c  "  sad  L=b  "  sdb  Lac  ~  |  gaoLtbd]'  \  SabLCcd]=  0  '  (l-130) 

Bab(Lde)  =  0  '  (1-131) 

In  chapter  2,  we  will  go  through  a  derivation  of  the  BMB  field  equations, 


from  a  slightly  different  Lagrangian  .  The  field  equations  here  can  be 

(21) 

found  by  work  analogous  to  that  in  chapter  2 

This  is  an  underdetermined  set  of  field  equations  (coordinate 

transformation  variations  produce  four  identites) ,  and  further,  the 

,  ,  Ta 

relation  between  g  ^  and  L^c 


is  more  complicated  than  the  one  in 


•  »•>  ■■  r;  : 


. 


general  relativity  ,  which  is: 

,d 


^ab;c  ^ab,c  ^ad  ^bc  ^db  ^ac  ^ 


(1.132) 

To  get  a  usable  set  of  equations  we  perform  a  "projective  transformation" 

rtc  =  +sbv0  d.133) 

where  V  is  any  vector.  Such  a  transformation  preserves  geodesics  and 
c 


parallel  transported  directions  globally 


(22) 


However,  the  lengths  of 


vectors  undergoing  parallel  transport  are  changed.  We  choose: 

vc=  3  L&*d]  (1,:L34) 

( .specifically  to  setf^^equal  to  zero).  With  this,  our  field  equations 


become: 

(1.135) 

(1.136) 

(1.137) 

(1.138) 

The  underdetermination  of  the  system  rests  in  eq.  (1.138),  in  which  V  is 
not  fully  determined.  Because  of  this,  (1.138)  is  usually  replaced  by 


^ab,c  ^ad  ^cb  ^db^ac  ^ 

r°  =  0 

[acj 

R/  ,  x  =  ° 

(ab) 

R  _  -  =  V  ,  -  V, 

[ab]  a,b  b,a 


the  equivalent  statement: 

RC(a.bLc]=  °- 

With  this  system  we  can  write  R  ^  as: 


Rab  =ftb,c  '  Kr(a0),8  +r(bc),a)  ^b  +  O' 


4 

cd 


(1.139) 


(1.14Q) 


3/ 

to  explicitly  display  the  symmetry  in  a,b  that  P (ac)  ^  now  has.  We  note 
that  (1.134)  is  equivalent  to  : 

(f£ifahl)  -  =  0  .  d-141) 


,6 


Eq.  (1.135)  is  not  quite  the  same  as  (1.132).  The  indices  in  the  second 


K. 
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term  connection  have  been  reversed,  so  that  in  fact: 
gab;c  2  H 


ad 1  Ccbl 


(1.142) 


If  we  want  the  covariant  dervative  of  g  ,  to  be  zero  we  must  set 

ab 


__  u. 

I  j~c^j  =  0.  Then  the  connection  is  symmetric,  there  is  always  one 
symmetric  g  ^  which  solves  (1.135)  (now  equivalent  to  (1.132)  ), 

and  we  find  the  field  equations  are  those  of  general  relativity  without 
electromagnetic  field.  ,  It  appears  that  g^^  is  necessarily  nonzero  in 
nonsymmetric  unified  field  theory,  although  this  will  lead  to  difficulties 
in  giving  the  metric  a  full  geometric  interpretation  (see  chapter  3) • 

We  have  not  yet  identified  the  representative  of  the  electromagnetic  j 
field  in  this  theory.  We  can  expect  some  modification  of  Maxwell's 
equations,  so  we  must  anticipate  a  generalized  electromagnetic  field 
which  approximates  EMT's  in  appropriate  circumstances.  Simply  on  the 
basis  of  numbers  of  components  we  could  say  that  g  is  either  the 
electromagnetic  field,  or,  on  the  same  grounds,  its  dual.  The 
field  equations  are  not  much  help,  unfortunately:  g^^has  the  divergence 
equation  (1.141 )  but  is  the  object  with  a  potential,  V^. 

In  the  absence  of  immediate  identification  of  the  electromagnetic 
field,  we  can  look  for  it  by  examining  the  properties  of  the  field 
equations.  In  particular,  we  can  look  for  solutions  representing  single 
isolated  charged  particles.  In  the  simplest  case  we  ask  for  a  spherically 
symmetric  point  mass  geometry,  asymptotically  Minkowskian: 

Kf~s(ab)  =  ?at  •  (1-143) 

For  the  elect ro magnetic  field,  we  seek  some  tensor  quantity  related  to 

a-  which  goes  nicely  to  zero  as  r  tends  to  infinity .  Let  us  not  be  as 

Cab} 

presumptuous  as  to  ask  for  a  regular  solution  yet.  This  type  of  solution 
would  be  comparable  to  the  Reissner- Nordstrom  geometry.  Unfortunately 


I 

' 

■ 

. 
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neither  the  use  of  gj-ab-j  nor  its  dual  as  the  electromagnetic  field  yields 
any  such  easily  interpretable  solution^”2^  .  This  fact  led  Wyman^2^ 
to  suggest  that  the  line  element  might  not  be  so  directly  related  to  g 

3j  D 

as  we  have  assumed.  This,  however,  is  an  intuitively  unappealing  thought. 

Another  approach  to  the  electromagnetic  field  is  to  take  the  weak 
field  limit  (ie,  linearization).  Since  EMT  is  valid  for  laboratory 
experiments  on  earth  (a. weak  field  situation) ,  we  expect  to  find  the  EMT 
equations  in  this  limit.  This  in  not  the  case:  the  linearized  equations 
governing  gf ,_are  weaker  than  Maxwell's  equations^).  Further, 


Papapetrou 


[ab] 

(23) 


showed  that  in  second  order  the  modification  of  R 


(ab) 


due  to  g  ^  is  not  in  any  way  similar  to  that  of  the  stress-energy - 
momentum  tensor. 

Yet  another  way  to  identify  the  electromagnetic  field  is  by  studying 
the  equations  of  motion  of  isolated  particles.  This,  however,  is  a 
difficult  problem  in  itself.  Being  a  large  and  complicated  subject, 
we  will  not  delve  into  the  mathematical  details  here.  In  general 
relativity,  the  one  body  problem  is  nontrivial,  and  the  two  body 
problem  is  tractable  only  in  approximation.  Collections  of  particles  in 
numbers  large  enough  to  be  approximated  as  fluids  are  treated  with  the 
stress-energy-momentum  tensor,  T  in  Einstein's  equations.  Individual 
particles  in  a  many-body  problem  are  viewed  as  localities  with  nonzero 
(mass -distribution-type)  T  with  vacuum  (T^  not  necessarily  zero) 
separating  the  individuals.  Unified  field  theories  prefer  not  to  have 
a  stress-energy-momentum  tensor,  so  equation  of  motion  work  is  usually 
done  with  isolated  singularities.  Most  researchers  rationalize  that 
there  should  exist  regular  solutions,  with  form  as  yet  unknown, 
asymptotically  equivalent  to  singular  solutions  (all  calculations 


' 


1 
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are  done  in  far  fields) .  Considering  the  scarcity  of  known  completely 

regular  vacuum  solutions,  the  author  believes  that  this  attitude  is  hardly 

justified.  It  is  possible  to  perform  this  kind  of  analysis  in  general 
(27) 

relativity  ,  but  apparently  the  state  of  the  art  of  equations  of  motion 

was  not  sufficiently  advanced  when  Weyl's  and  five-dimensional  theories 

were  being  actively  studied,  as  very  little  such  work  was  done  in  them. 

Early  work  in  the  nonsymmetric  theory suggested  that  it  did  not 

have  the  expected^  ^  equations  of  motion,  but  recent  work^^  has  been 

(31) 

more  encouraging.  Johnson,  in  a  series  of  papers ,  develops  a  very 
general  approximation  applicable  to  the  study  of  equations  of  motion. 

From  the  application  of  his  approximation  to  Einstein's  theory  he 


•  cL*b 

identifies  the  electromagnetic  field  as  the  dual  of  g  : 


^ab  2eabcd  ^ 


.(pd] 


(1.143) 


where  e  is  Levi-Civita's  tensor  density.  He  states  that  this  is  the 
simplest  expression  approximating  the  field  which  governs  the  motion  of 
charged  particles  at  intermediate  (laboratory)  distances,  whereas  the 

(32) 

true  field  governing  motion  everywhere  is  of  complicated  form  .  With 
this  identification,  simple  charged  "particles"  interact  via  the  Lorentz 
equation  of  motion  with  radiation  reaction  over  intermediate  distances. 

The  interaction  at  microscopic  distances  is  complicated  and  not  yet 
resolved;  the  interaction  at  astronomical  distances  involves  an  asymptotic 
field  that  does  not  vanish  at  infinity . 

This  appears  to  be  the  state  of  the  standard  nonsymmetric  theory  today. 
Work  is  still  being  done,  particularly  by  Johnson,  who  is  still  applying 
his  approximation  scheme  to  the  theory.  One  variant  of  the  nonsymmetric 

.  (33) 


theory,  other  than  the  BMB  version,  is  Kursonoglu  s 

(34) 

summarize  his  theory.  He  represents  s 


We  will  just 


-  < 


. 
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^ab  ®(ab)  +  q  ®£ab] 


where  g,  x  is  interpreted  as  the  metric  tensor,  g 


(1.144.) 

as  the  generalized 


’(ab)  ~  & [ab] 

electromagnetic  field,  and  q  is  a  constant  with  the  dimensions  of  an 

electric  field.  From  a  variational  principle  he  derives  field  equations 

which  can  be  put  into  the  form: 

R 


(ab)  =  ^ 


-det  g 


-det  g 


(g(ab)  -  ‘  SCaU 


ab 


R[[ab],c]  +  2 ^  S[[abJ,cJ  0 


£ 


Cab} 


=  0 


pd  pi  d 

^ab , c  ^ad  '  cb  ^db '  ac 
where  K  is  related  to  q  by: 


0 


(1.145) 

(1.146) 

(1.147) 

(1.148) 


K 


4G 

4 


(1.149) 


[ab] 


is  the  antisymmetric  part ,  but 


cL*b 

Here  g  is  the  inverse  of  g  ,  and  g 

3/  0  | 

(ab) 

7T  )  is  the  inverse  of  g^.^  .  Eq. (1.147)  is  again  equivalent  to  de¬ 
manding  that  torsion  be  traceless.  Kursonoglu  claims  to  have  correct 
equations  of  motion  and  a  class  of  completely  regular  solutions  which 
have  complicated  magnetically  charged  structure. 


Chapter  II 

Modification  of  the  Nonsymmetric  Theory 
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History 

If  we  are  discouraged  by  the  behavior  of  Einstein's  nonsynrmetric 
theory  (even  today  the  best  conceivable  case  implies  that  charged 
particles  have  a  nonzero  asymptotic  electromagnetic  field) ,  we  have  the 
option  of  modifying  it  in  the  hope  of  finding  better  properties.  In  1954, 

.(35) 


Bonnor 


proposed  modifying  the  Lagrangian  in  Einstein's  variational 


principle  to  improve  the  equations  of  motion.  His  variational  principle 


was  equivalent  to: 


$/(,/=? gaV.  +  ^sCa-bl 


ab  T  2 
k 


Lab] 


%  4 

)  d  x  =  0 


(2.1) 


where  k  is  a  universal  constant  with  dimensions  of  length.  With  this  it 
was  found  that  the  equations  of  motion  became  acceptable  in  a  slow 
motion  approximation  scheme.  He  proposed  the  interpretation  that  the 

P  a  "hi 

electromagnetic  field  is  the  dual  of  g  (as  G.R.  Johnson  does  in  the 
usual  theory  (eq  (l .  142))) ,  but  was  unable  to"  find  a  reasonable  point 


charge  solution. 


(36) 


Moffat  and  Boalw~'  (1974)  took  the  interpretation  that  g  ^  itself 
is  proportional  to  the  electromagnetic  field  and  used  the  Lagrangian  of 
eq.  (2.1)  to  find  the  field  equations  (we  will  go  through  the  derivation 
shortly) : 

s  -g  Pd_gnnd=o  (2.2) 

gab,c  gad  '  cb  gdb  ac 

(2.3) 

(2.4) 

(2.5) 


r, 


[ac] 


=  0 


R(ab)  +  1  (ab)  =  0 


R[ab]  +  1  [ab]  L , b  '  Vb,a 


_  -4tt 


'ab 


2  ^Sac 


.CcdJ  i  Ccd) 

Sdb  2  gab  g[cd]g 


+  g[ab]  ^ * 


(2.6) 


where : 


28 


Compare  this  set  to  eq.  (1.135)-(1.138) ,  in  the  same  notation's  the 
modification  is  the  introduction  of  I  .  Again,  we  prefer: 

ci  0 

[(RCabj  +  1  [ah]  \cj  =  0 

to  eq.(2.5),  since  the  potential  V  is  not  determined  by  the  field 

Si 

equations.  Here  eq*,,  (2.3)  is  equivalent  to:. 

(  /TgIaM)  >b  =  0  .  (2.8) 

Note  that  because  the  Lagrangian  now  has  two  terms,  the  second  being 

electromagnetic  and  having  a  new  constant  k,  this  theory  is  not  as  nicely 

unifiedas- 'Einstein's  nonsymmetric  theory.  We  hope  that  this  will  be 

offset  by  better  behavior  of  the  new  theory. 

Although  Bonnor's  treatment  of  the  equations  of  motion  also  applies 

(37) 

to  this  identification  of  F  in  a  more  modern  analysis  Johnson 
showed  that  in  his  approximation  scheme,  charged  particles  interact  in 
first  nontrivial  order  according  to  Maxwell's  equations  and  the  Lorentz 
force  law  with  radiation  reaction,  but  without  the  asymptotically 
nonzero  field  found  in  Einstein's  UFT. 

Moffat  and  Boal  found  a  static,  spherically  symmetric,  electrically 
charged  solution,  which  looks  like  a  generalization  of  EMT's  Reissner- 


^(38), 


Nordstrom  solution 


1  -  2m  + 


47^ 


dt2  - 


dr 


(l  -  2m  +  4ttQ  ) 


,  ,  -  7  (2.9) 

-  r  (  d©^  +  sin  9  d^  ) 

where  m  is  mass  and  Q  is  charge.  The  only  nonzero  element  of  g|-ab-j  in 
the  solution  is: 

e  =  M.  (2.10) 

gri41  r2  • 

On  the  basis  of  this  we  take  the  interpretation  that  g^^is  the  electro¬ 


magnetic  field  multiplied  by  k: 


' 

, 
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SCaM=kFab  '  (2-11) 

Then,  if  k  =  0,  this  solution  becomes  exactly  the  Reissner- Nordstrom 
solution,  suggesting  the  possibility  that  the  field  equations  might 
become  EMT's  in  the  limit  where  k  tends  to  zero.  In  fact  this  is  the  case 
if  we  take  the  limit  of  (2.2) , (2.4) , (2.5) , (2.8) : 

=  0  (2.12) 

(2.13) 


3ab;  c 


R  ,  =  8ttT  , 
ab  ab 


-8ir  F  =  V  ,  -  V, 

ab  a , b  b , a 

(/=I'Fab)ilj  =  0 


Tab  ^a^db  +  ‘tSabFcdF 


cd 


(2.14) 

(2.15) 

(2.16) 

Taking  the  limit  is  not  quite  a  straightforward  procedure.  The  small  k 

limit  of  eq.(2.8)  is  a  trivial  equation,  because  there  is  an  overall 

factor  of  k.  In  order  to  get  the  electromagnetic  divergence  equation  we 

must  ask  for  a  nontrivial  equation  in  the  limit  (one  way  is  to  say  that 

k  is  very  small  but  not  zero).  Nonetheless,  this  limit  is  an  interesting 

property  not  found  in  Einstein's  UFT.  From  eq.(2.14)  we  identify  the 

"generalized  electromagnetic  potential"  vector,  A  ,  by: 

a  i 


V  =  -8ttA  , 

a  T  3. 
k 


(2.17) 

although  in  the  k  ^  0  case  it  is  not  directly  related  to  g  >  tut 
rather  to  R 

The  line  element  stated  in  eq.(2.9)  is  in  the  real  metric 
interpretation.  We  can  do  the  same  work  with  a  hermitian  symmetric 
metric  obtaining  the  results  that  k  becomes  imaginary  (ik)  and  the 


TabJ  ‘ 


solution  line  element  becomes: 

is2  =  (1  -  (1  -  2m  +  ^f-)  dt2  -  - -  2 

r  r  r  (l  -  2m  +  4~ffQ  ) 


dr 


. 


. 
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2  /,  2  .  2_  , /2\ 

-  r  (d©  +  sin  ©  dp  )  . 


(2.18) 


Inside  the  radius  kQ  the  signature  of  the  metric  changes  from  the  usual 

(+ - )  to  ( - ),  which  is  a  pure  spacelike  signature,  so  no  particle 

following  a  timelike'  path  can  ever  penetrate  into  this  region,  in 


particular  to  get  to  the  central  singularity.  Moffat 


(39) 


has  succeeded 


in  showing  that  this  solution  is  world- line  complete,  and  since  the 

singularity  is  not  accessible  to  real  particles,  the  physical  space 

is  completly  regular.  Particles  which  penetrate  to  radius  kQ  are  not 

deflected  back  into  the  same  space,  but  move  onto  a  second  Riemann 

(40) 

sheet  of  the  maximally  extended  geometry  '  .  Having  such  a  "regular" 
solution,  the  hermitian  symmetric  alternative  may  prove  to  have  some 
advantage  over  the  real  metric  theory.  We  will  retain  the  real  metric 
interpretation  (unless  otherwise  stated)  for  reasons  discussed  in 
chapter  3* 

For  discussion  of  magnetic  monopoles  see  "Duality  and  Magnetic 
Monopoles",  later  in  this  chapter.  In  other  work  done  to  date  the  first 
steps  towards  finding  a  charged  spinning  (analogue  of  the  Kerr-Newman  ) 
solution  have  been  discussed^"^  ,  and  several  extensions  of  the  theory 
have  been  proposed . 


Derivation  of  the  Field  Equations 

Let  us  now  derive  the  field  equations  from  a  variational  prin¬ 
ciple^2'*'^.  We  choose  the  Lagrangian  in  analogy  with  EMT  as  follows: 


y?sa\b  - 


[ab]  , 


(2.19) 


where  the  contracted  curvature  is  denoted  by  B^,  and  the  connection 

ebb  Si 

by  as  in  eq.(l .129)-(l .131) •  The  quantities  g  and  are  to 
be  varied  independently.  We  have  written  the  formal  variation  several 
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times;  taking  the S  inside  the  integral: 

/kb  s(/V  £ab  ^  -  Wabl 

-y  L  k 


S  Cab]  ^ 


4 

d  x  =  0. 


Let : 


-  H 


Cab])  =  Lb 


so  that : 


kb +  Lb)5 


ab  ab  0  .4  _ 

i  +  £  °  Bab  d  x  =  0 . 


Since  depends  only  on  the  connection,  its  variation  is 


of  £gab  ,  while  Sgab  is  obviously  independent  of  &L^  : 


ab  . 


(Bab  +  Lb  )  5  ^  ^  =  0 


,4 


rab  Sb  ,  d'x  =  0 


ab 


Eq.  (2.23)  implies: 

B  ,  +  I  ,  =  0 . 
ab  ab 

Let  us  now  determine  I 


ab  ' 


S(-Srab]S[abJ)  =  8[WJSfiab  +  £rab4« 


’ba 


Now: 


5sad  =  'Sbd  Sac  °S 


C  be 
0{ 


and  from  properties  of  determinants: 

£  £  ab 

ag  =  -ggabdS 


so : 


C  ab  , — ,c  db  1  ab  C 

dg  =  /-gdg  +  -  g  og 

2Vcg 

r — 1 £  ab  C  ab  1  / — >  ab  c  cd 

f-s  dg  =  dg  +  i Ag  g  gcd  £  s 


Then: 


)  =  Sr.JsrabJ  -  Sar,  s,v 


[ab] 

^ fab] J  e£ba]va  &ac  °db 


ab 


£ 


4/^Sd/°dJgab  S, 


[ab] 


(2.20) 

(2.21) 

(2.22) 

independent 

(2.23) 

(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 


(2.31) 


Now: 


. 


1 


(2.32) 


fs  gab  Se 


ab 


=  — cr 

o 


ab 


s  ga6) 


implies: 


Cab]  \  C  ab  fed]  C  ab 

*r-&  Sr^i)  =  grba3«fi  -  Sa„g  SrlA 


Tab] 

i  Ccd]  C  ab 

+  *  S  ScdSatM 


ac 


’db 


(2.33) 


Therefore: 


'ab 


-4~r 

k2L 


[cdj  i 

SacS  sdb  2  SabS[cd]S 


Ccd] 


+  g 


Cab] 


(2.34) 


For  the  variation  of  the  curvature  we  use  the  variation  of  the  co¬ 
variant  derivative  of  the  connection: 

-  L*,  &l£ 

cd  De 

to  find: 


<s^,d  =  Sl?=,d  +  Ltd  K  - 


§B  ,  =  -$LCV  +  ,  +  L~  dL~  +  IT,  dV 

ab  ab,c  ac,b  ac  db  db  ac 


■d  §tC 


c  §Td 


bdc  bab 


Slid,  -  Ld  Si, 


d  §L" 

ab  dc 


=  -(SltV  +  (5l°)  :  2L 


ab  ;  c 


ac'  ;b 


d  £tc 
febj  ad’  * 


Using: 


La  LCac] 


0  =  fgab  (-(  Sl°J  +  (  <§LC  )  +  2L  d  SlC  )  d4x 

J  a  v  ab2 ;  c  v  ac';b  £cbl  ad' 

we  find,  after  some  algebra: 

b  cb  /  ad 

Jfcdj  c  £  ;d 

,4 


0  ab  ,  ,  _  ad  T  b 

-  2  g  L  +2g  L, 

J  C  ^  C 


&  L°  d  x  . 
ab 


Define: 


^ab  ab 
P  c  =  S 


b 

[cd] 


(2.35) 


(2.36) 


(2.3?) 

(2.38) 


(2.39) 


(2.40) 


i  0  ab  2  Cb  ad  9  ad 
,c  -  %  Lc  +  3^  Ld  +  2£  L 

and  (2.39)  becomes: 

>abo  -  s>add)  K/x  -  0  . 

Therefore  pab  must  vanish.  With  a  little  calculation  we  can  find  the 
c 

equivalent  statement  in  g^  c  rather  than  g  ^ c : 


(2.41) 


. 
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0  gab,c  gadLcb  SdbLac  3  gac  L[bdl  3  SabLCcdl  (2.42) 

which  we  have  already  seen  as  eq.(l.l30). 

Now,  as  discussed  in  chapter  1,  we  perform  a  projective  transformation: 
a  _a 


rbc  =  +  f  SbLc 


(2.43) 


(which  makes  the  new  torsion,  (V-^j ,  tracelesss).  It  is  fairly  easy  to 


see  that: 


R 


ab 


=  B  -  f  (L  ,  -  L,  ) 
ab  3  a,b  D,a 


and  that  (2.43)  and  (2.42)  respectively  are  equivalent  to: 


ab 


db  na  ,  aanD  ±  ao  ,  pa  n a  \  A 

a  i0  +  s  rdc  +  s  rcd  - 1  £  ( rdc  +rod),_  =  o 


ab 


(2.44) 

(2.45) 

...  ,  0  ,  ,  .  -  .  (2.46) 

db  '  ac  ad  1  cb 

Gathering  everything  together  our  field  equations  are  eq.(2.46)  together 


gab,C  grf-h  P,  ’  ® 


with: 


R(ab)  +  X(ab)  ° 

RCab]  +  Jrabj  “  3  ^a,b  "  Va^ 

r  c  =0 
1  CacJ 


(2.47) 

(2.48) 

(2.49) 


where  I  T  is  determined  by  (2.34).  Now,  L  is  not  fully  determined,  so 
ab  a 


we  merely  identify  a  potential  vector  for  R 


£abj' 


V  =  %  L 

a  3  a 

RCabJ  +  LCabJ  =  Va,b  '  \a  ' 


(2.50) 

(2.51) 


This  completes  the  derivation  of  the  basic  field  equations,  (2.7) 
following  immediately  from  (2.51).  To  verify  (2.8): 

(£&0:l))C  =  0 

we  subtract  the  contraction  of  (2.4 5)  on  a  and  c  from  the  contraction  on 


b  and  c: 


0  =  2  £ 


facj 


=  2  £ 


fac3 


,c 

+  2  £ 


1  a-3'6-  r  ° 

2  a  1  cd 


1  da  nc  _  i  p 

2  £  /  2  £  I 


i  __ad  nC 


dc 


dc 


+  f  ^  r° 


cd 


(ad)  n  c 
cd 


(2.52) 
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o  =  /ac3 


»c 


(2.53) 


Note  that  covariant  differentiation  with  repect  to  L^c  was  used 
significantly  in  this  derivation  (in  steps  done  to  find  eq.  (2.39)). 

Oo variant  differentiation  with  respect  to  is  not  the  same  operation. 

Physical  interpretation  will  require  that  we  single  out  one  covariant 
derivative,  and  the  choice  made  may  he  hard  to  justify.  Einstein  and 
Straus  were  able  to  derive  their  field  equations  (1 .135) -(1 .138) 

from  a  Lagrangian  without  a  projective  transformation,  by  including 


extra  terms  in  the  Lagrangian  with  Lagrange  multipliers  to  obtain  the 

r„h-T 

traceless  torsion  and  divergence  equation  on  g  .  With  a  similar 
derivation  in  BMB  theory  we  would  not  need  to  consider  L^  .  We  will 
return  to  consideration  of  covariant  differentiation  in  chapter  3* 
Duality  and  Magnetic  Monopoles 

It  is  wall  known  that  Maxwell’s  vacuum  equations  are  dual  invariant, 
that  is,  if  F  ,  is  a  particular  solution  of  them,  then  another  solution 

cl  D 


is: 


F  ,  =  F  ,  cos«+  F  ,  sin  <=* 

ab  ab  ab 


(2.54) 


* 


where  o<  is  an  arbitrary  angle,  and  F  is  the  dual  of  Fab  : 

C^abcd^ 

(  £  is  the  Levi-Civita  alternating  symbol).  Eq.  (2.54)  is  called  a 
abed 

"dual  rotation"  because  it  is  reversible  in  the  sense  of  a  rotation: 


i  '  * 

F  ,  =  F  ,  cos  (-<*)  +  F  ,  sin  (-<*). 
ab  ab  v  '  ab 

The  key  relation  that  makes  eq(2.56)  work  is: 

**  _  i  £  F*cd 
ab  2  ^abed 

=  ic  <7  ?7Ce  „df  f 

*  ^abed  efhj  u  7?  (  '  kl 


(2.56) 


=  SCa  SbJ  Fkl  =  'Fab  • 


(2.57) 


The  dual  of  the  dual  is  the  negative  of  the  original  field.  If  we 
relax  the  usual  source  equations  to  allow  the  possibility  of  magnetic 
charge,  the  complete  Maxwell  equations  (with  sources)  become  dual 
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invariant,  and  such  rotations  (by  o(=  )  generate  magnetic  monopoles 

out  of  electric  monopoles.  Dirac  showed  that  a  viable  electromagnetic 
theory  can  be  constructed  with  such  monopoles,  leading  to  charge 
quantization^^  . 

EMT's  (vacuum)  equations  are  similarly  dual  invariant,  with: 


ab 


“  If 

e  ,  ,  F 

abed 

1 

ce 

2 

eabcd  S 

(2.58) 


ef 


where  e  ,  ,  is  the  alternating  tensor: 

abed 


eabcd  ^  fabed 


(2.59) 

(2.60) 
(2.61) 

Therefore  every  electrically  charged  monopole  solution  has  a  corresponding 


4  eabcd  eefhj 


eabcd;e  ^ 


ce  df  hk  jl  rl  rk 
g  g  g  g  = 


Dirac  monopole  solution.  The  general  electrically  and  magnetically 

(44) 

charged  Reissner-Nordstrom  solution  is  thenv  : 

.2  _2 


ds2  =  (1 . 25  +  2+  L-)j  at2 

-  r2(  d92  +  Sin2©  d{52  ) 


(2.62) 


dr 

2  2 
(l  -  2m  +  4tt(Q,  +  L  )) 

r  2 

r 


0. 


E  =  9 

r  2 
r 


B  = 
r  2 
r 


(36) 


(2.63) 

In  proposing  their  theory  Moffat  and  Boalww  suggested  that  their 
equations,  having  an  electric  monopole  solution,  probably  do  not  have  an 
equivalent  magnetic  monopole.  If  this  were  indeed  the  case,  it  would  be 
a  favourable  property  of  the  theory ,  as  no  magnetic  monopoles  have  yet 
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been  detected  in  nature 


(^3) 


No  such  solution  has  been  found  in  BMB 


theory  by  studying  plausible  trial  solutions  involving  polynomial 

(45) 

expansions  in  the  radial  coordinate  ' , 

Let  us  consider  dual  rotations  in  BMB  theory.  The  lack  of  an 
obvious  fundamental  tensor  for  raising  and  lowering  indices  (see  chapter 
3),  the  inclusion  of  the  electromagnetic  field  in  the  metric,  and  the  lack 
of  any  covariantly  constant  tensor  related  to  the  metric  make  it 
difficult  to  define  a  usable  duality  rotation.  In  analogy  with  EMT  we 
formally  write: 

s[ab]  =  SCab:  00S*  +  <£abj  sin«'  (2'64) 

There  are  a  number  of  fairly  simple  possible  forms  to  choose  from  for 
* 

the  dual,g 


[ab]  ‘ 


g 


3 


g 


* 

_  i 

_[cd] 

[ab] 

2  abed 

g 

* 

„  1 

ec  df 

[ab] 

2  eabcd 

s  g  g[ef] 

* 

__  1 

ce  df 

[ab] 

2  eabcd 

g  g  gref] 

* 

(ce)  (df ) 

[ab] 

—  —  r 

2  abed 

g  g  g 

* 

__  1 

ce  df 

[ab] 

2  eabcd 

r  ®  [ef  ] 

[ef] 


g 
ab  . 


where  TT  is  the  symmetric  tensor  satisfying: 

^ab  ra 

*  s(cb)  _i>0  ’ 


and: 


6  * 

s[ab] 


V  e 


2  abed 


4 


cd 


where : 


7ab 
<t>  g 


- 

[cb]  “  ^c 


(2.65) 

(2.66) 

(2.67) 

(2.68) 

(2.69) 

(2.70) 

(2.71) 

(2.72) 


(fh  is  antisymmetric) .  If  we  decide  gab  is  the  fundamental  geometric 

1  "X-  C^d3 

tensor(see  chapter  3)  »  then  g^^is  the  most  natural  choice,  but  g 

does  not  formally  look  like  g with  two  indices  raised.  Eq4(2.66) 


w- 


■ 
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(2.69)  seem  closer  formally  to  the  EMT  dual,  ^g  j-^-j  being  the  most 

4  * 

likely  if  ls  the  fundamental  tensor  (and  g  possible  in  this 

2  * 

case).  Eq.(2.66)  and  (2.57)  are  alternates  to  (2.65).  g  j-^-j  has  a  more 

cl/b 

natural  order  of  summation  indices  in  view  of  definition  of  g  ,  but  the 

"upper  index  field"  is  not  antisymmetric.  In  ,  the  upper  index 

6  * 

field  is  still  skew.  g  might  be  chosen  if  nothing  else  is  suitable. 

Let  us  look  in  more  detail  at  each  of  these  proposed  duals .  Under 

6  * 

scrutiny  g  j-^-j  can  quickly  be  eliminated.  The  determinant  of  g can 
be  zero  (example:  the  electric  monople  exact  solution),  so  ^g  does 

not  always  exist.  Further,  when  it  does  exist,  we  find: 

*  ~  d6t  S  [at] 

/ab  -  EE  .  abed 


g 


*  zip 
_  /S 

Cab)  Cab] 


[cd] 


(2.73) 

(2.74) 

(2.75) 


This  does  not  become  the  usual  dual  in  the  small  k  limit.  To  investigate 


1  * 
g 


Cab]’ 


[ab]  ,  ...  (46) 

can  be  written 


[ab]_  -ab  +  ^ac^td 
-g  r  -g 


[cd] 


(2.76) 


(2.77) 


where : 

**  S  det  s(at) 

when  </>  is  nonzero  (^  zero  must  be  treated  separately) .  With  this: 

)  ,  (2.78) 

giving  a  dually  transformed  field  (replacing  the  prime  with  the  number 

of  the  dual  chosen) : 


1  *  _  (5* 

S  [ab]  g  [ab]  -g  1  S  [ab] 


1 


S[ab]  S[ab] 


+  ^  <*. 


(cos«  +  ^  si n°< ) 

)  sin  °< . 


ce  _df 

abed  ^  *  SCef] 


(2.79) 


This  implies: 

V  -  det  (  g  £abj  ) 


(2.80) 


. 
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T6  ^ 
1 


g  -  det  g 


ab 


(2.81) 


/  g  • 


The  supposed  inverse  transformation  is  complicated: 

)*  = 


1,1 

(  g 


[ab] 


*  _/jT  A  \  M  V  ,  i 

" /r 7  ( scabi}  +  r  (  * 


'abed 


ce  df  /I  \ \ 
r  r  ^  g  [ef  J 


(2.82) 


Q 


Cab] 

=  gEab] 


g[ab]COS<*  (  g  Tab]  ^  sinc?< 

2 . .  ^  _  yy 


cos*c<  +  (*£%  -  )  coso<  sinc<  sin2c* 

^  /ig 


5  * 
+  g 


tab +  W=P  sin2*+rcos*  sin*  A  -  j-) 
L  -g*-g  -g/-g  s  -g  ^ 


(2.83) 


5  ** 

-  g 


[ab]  1 


2^  2 
6  sin  « 


gg 


where : 


5  **  =  i  e  e  ^e^f  kh  ij 
g[ab]  4  eabcd  efkl  ^  V  ^  ^  s£h,j]  * 


(2.84) 


We  have  certainly  not  been  returned  to  the  original  field  and  therefore 
1  * 

g  j-abj  does  not  generate  a  dual  rotation.  (Since  the  <j>  nonzero  general 
case  fails,  we  will  not  worry  about  the  special  case  of  zero  <f) . 

Note  eq.(2.84):  it  is  a  "dual  of  the  dual"  expression  analogous  to 
eq.(2.56).  This  expression  (in  case  (4)),  or  a  similar  expression  >  appears 
in  the  inverse  transformation  equations  for  all  of  cases  (2)  through  (5) . 
One  of  the  conditions  neeessary  to  obtain  a  dual  rotation  is  that  the 
analogue  of  (2.60)  (which  implies  (2.57))  hold,  and  this  appears 
difficult  to  satisfy.  In  fact,  there  are  other  nontrivial  problems 
in  cases  (2) ,  (3)  and  (4) ,  but  this  is  the  only  bistacTe'  in  case  (5) . 


It  is  necessary  that ; 

abed  n 

y-.  =  U 


(2.85) 


.(47)' 


(for  example,  this  defines  how  tensor  densities  transform^  ) ,  and  the 
upper  and  lower  alternating  densities  are  related  by  their  identities, 
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including: 

i-£abcd  e 


_  r  a  rb 

cdef  “  d[f  de]  ' 


(2.86) 


Taking  the  covariant  derivative.. of  eq.(2.86)  we  find: 

^abcd;e  ^  (2.87) 

Now,  because  the  covariant  derivative  of  (and  hence  ?rat)  is  not 

zero ,  obviously : 

-WVJ  y  s*fS ^  (2.88) 

and  we  conclude  that  none  of  the  possibilities  we  have  suggested  generate 
dual  rotations . 


**  6abcd  eefkl 


From  the  considerations  of  this  section,  it  appears  that  there  is  no 
simple  way  to  define  a  meaningful  electromagnetic  dual  in  nonsymmetric 
unified  field  theories,  much  less  test  dual  invariance.  While  experiment 
presently  indicates  no  symmetry  between  electric  and  magnetic  charges, 
the  acceptability  of  the  total  Iobs1  of  the  duality  (.concept  is  a  matter 
for'  investigation,’ 

Possible  Special  Case  Solutions 

This  section  points  out  some  special  circumstances  in  BMB  theory  in 
which  solutions  may  exist .  These  cases  are  special  because  they  may 
shed  more  light  on  the  differences  between  BMB  and  EMT,  and  possibly  on 
the  validity  of  BMB  theoiy  itself.  Assumptions  of  a  few  simple  particular 
forms  has  not  led  to  any  such  solutions  yet. 

An  important  question  to  ask  is:  are  there  any  solutions  involving 
nonzero  electromagnetic  field: 

®Cab]  *  0  (2'89) 

on  a  flat  background?  We  will  have  to  specify  what  we  mean  by  "flat". 

We  can  ask  for  a  Minkowskian 
®(ab)  ^ab 


(2.90) 


' 


or  that  there  exist  coordinates  in  which  geodesics  are  straight  lines: 


r, 


=  0 


(be) 

(everywhere)  or  that  there  be  no  curvature: 
^bed  =  0  * 


(2.91) 


(2.92) 

These  are  non- equivalent  possibilities,  becoming  equivalent  in  the  small 
k  limit.  Eq,  (2.92)  is  certainly  the  strongest  condition.  (  No  such  flat 
background  solutions  exist  in  EMT:  nonzero  Maxwell  field  implies  nonzero 
stress-energy-momentum  tensor,  and  hence  nonzero  Ricci  curvature.  The 
BMB  equations  become  EMT's  as  k  tends  to  zero,  so  we  can  immediately 


say  that  any  solution  with: 

tab]  ~) 

*  / 


lim 

k-K) 


1  0 


(2.93) 


violates  all  of  (2.90)  -  (2.93)  because  of  the  small  k  limit. 

This  last  statement  raises  a  more  general  question:  are  there  solutions 
which  do  not  have  EMT  analogies?  Such  solutions  would  become  purely 
gravitational  as  k  tended  to  zero ,  having  all  electromagnetic  quantities 
functions  of  k  which  go  to  zero  in  the  limit .  If  there  is  the  possibility 
of  such  solutions,  then  I  ^  must  be  studied  to  see  if  there  are  any  circ¬ 
umstances  in  which  we  can  have  (2.89)  and: 


I 


(2.9*0 


I[Cab],o] 


(2.95) 


Only  in  such  cases  could  (2.92)  hold. 

Wyman  and  Zassenhaus^^  have  found  all  solutions  with  zero 
curvature  in  Einstein's  nonsymmetric  theory,  most  of  them  physically  non¬ 
trivial.  The  analogous  set  in  BMB  promises  to  be  reduced  by  the  constraint 


(2.94), (2.95). 

Another  (possibly  empty)  family  of  solutions  without  EMT  analogue 


would  be: 
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g 


Cab] 


(2.96) 


V  =  6 
a  ,a 

( j  is  any  sca.lar  field),  which  means: 


(2.97) 


R[ab]  _I  [ab]  *  (2.98) 

The  small  k  limit  implies: 

.  (2.99) 

k-X) 

Since  there  does  not  seem  to  be  any  obvious  interpretation  of  this 


family,  we  cannot  declare  its  members  desirable  or  undesirable 


without  further  investigation. 


Chapter  III 

Physical  Interpretation 

of  Geometry 
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General  Relativity's  Geometric  Structure 

General  relativity's  geometric  formulation  is  aesthetically  appealing 
because  almost  all  of  the  geometric  structure  has  physical  interpretation. 
There  is  very  little  extraneous  geometry,  v^ery  little  of  classical 
physics  left  out  (although  quantization  is  still  a  stumbling  block) . 

The  geometric  structure  is  a  four  dimensional  manifold  with  metric. 

The  four  dimensions  are  space  (three)  and  time  (one)  which  are  bound 
together  in  the  manifold,  spacetime.  The  metric,  g  ,  is  a  smooth 

ci  D 

(except  in  special  circumstances)  tensorfield  on  the  manifold  used  to 
obtain  the  "inner-product"  of  two  vectors: 


(3.1) 


and  the  "length",  a  special  case  of  this: 

lk!2  =  Sah  ka  •  (3-2) 

Now,  as  the  inner  product  is  usually  formulated  (based  on  its  beginnings 
in  ordinary  vector  analysis) ,  it  is  independent  of  the  order  of  its 
arguments,  and  length  is  always  non-negative.  Then  we  can  use  it  to 
define  the  "angle"  between  two  vectors  *  General  relativity  retains  the 
order  independent  nature  (this  means  g  ^  is  symmetric) ,  but  to  model 
light  propagation  properly,  the  non-negative  length  requirement  is 
relaxed.  To  have  length  of  consistent  sign  the  metric  must  have  consist¬ 
ent  signature:  (  ++++), or-  ( - ),  in  the  Sylvester  canonical  form^°\ 

Relativity  distinguishes  time  from  space  by  switching  one  of  the  four 
diagonal  elements:  we  choose  the  convention  (+---).  Then  there  are 
three  types  of  vector:  spacelike  ( |k|  negative),  pointing  "in  three  space" 
null  (|k|2zero),  pointing  along  paths  followed  by  light  in  four-space; 

O 

and  timelike  (|k|  positive),  pointing  along  massive  particle  paths. 

We  can  still  use  the  inner  product  to  identify  the  angle  between  two 


. 


spacelike  or  two  timelike  vectors,  for  example,  and  establish  other 
relations  for  other  cases.  Physically,  we  interpret  timelike  magnitude 
as  the  rate  at  which  time  passes  (as  already  mentioned)  for  an  observer 
on  the  curve  to  which  the  timelike  vector  is  tangent . 

On  a  manifold  an  affine  connection,  fbc  defines  quantities  such  as 
parallel  transport  and  covariant  derivatives .  The  covariant  derivative 


of  tensor  Q' 


3/  •  •  . 


b .  .  . : 


Q 


3/  ■  •  • 


b  .  .  . ;  c 

is  itself  a  tensor.  The  ordinary  derivative: 


Q' 


,3,  •  •  • 


b...,°  N  c 


w 


.3 •  %  •  • 


Id  •  •  • 


) 


(3-3) 


(3.4) 


is  not.  The  form  of  the  covariant  derivative  of  a  vector,  for  example, 

is: 

*VkVk°r:t  •  (3-5) 

Note  the  order  of  indices  in  the  second  term:  the  connection  is  not 
a  priori  symmetric.  To  parallel  transport  a  tensor  Q,  ’ ’ ’  along 

U  •  •  • 

3 I 

a  curve  with  tangent  k  ,  the  tensor  is  changed  as  it  moves  so  that  its 


covariant  derivative  in  direction  of  k  is  zero: 


Q‘ 


3/  •  •  • 


k  =  0  . 


(3.6) 


b  .  .  . ;  c 

Mathematically,  parallel  transport  is  thought  of  as  the  way  to  move  a 
tensor  taking  into  account  changes  in  local  curvature  without  changing 
what  the  tensor  actually  represents.  If  we  ask  for  those  curves  along 


which  the  tangent  is  itself  parallel  transported: 

ta  tC  =  0 
;c 

we  find  geodesics,  the  generalizations  of  straight  lines  on  curved 


(3.7) 


manifolds. 

An  important  theorem  is:  if  there  is  a  smooth  symmetric  two-index 
tensorfield  on  a  manifold,  in  particular,  the  metric,  gab,  then  there  is 


c 


■  :  -  r* 
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a  symmetric  connection  on  the  manifold  uniquely  determined  by  the  con¬ 
dition  that  be  preserved  in  parallel  transport,  ie: 

0  ~  Sab;C  d  d 

=  gab,c  "  sad  ^bc  "  gdb  ^ae 

(  an  equation  we  have  already  “seen) .  It  is  called  the  Levi-Civita  connec¬ 
tion,  and  has  the  form: 

^bc  =  i  s  (gdc,b  +  gbd,c  "  sbc,d^  *  ^.9) 

With  this  connection  inner  products  (and  hence,  lengths)  are  preserved  in 

parallel  propagation.  Inner  product  preservation  is  consistent  with  the 

intuitive  idea  that  parallel  transport  moves  vectors  maintaining  their 

intrinsic  nature.  Further,  it  is  easy  to  show  that  the  curves  that 

extremize  arc  length  between  pairs  of  points,  that  is,  which  are  mathema- 

(51) 

tically  the  solutions  ofVJ  ': 


ds  =  0 


,2  ,  a,  b 

ds  =  g  ,  dx  dx 
ab 


(3.10) 

(3.H) 

are  geodesics  of  the  Levi-Civita  connection.  This  reinforces  the  concept 
that  geodesics  are  the  generalizations  of  straight  lines. 

The  connection  defines  the  curvature  tensor: 


a  =  na  na 

bed  1  bd,c  bc,d 

and  its  contractions: 

Rab  ~  R  acb 

R  =  gab  R 


,  p3>  P®  _  Pa  p' 

'  ec  be  ed  be 


ab  ’ 


The  curvature  satisfies  several  identities,  in  general  relativity: 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

b  [cd;  e^ 

(and  others  when  a  is  lowered).  (3-17)  is  the  so  called  Bianchi  identity. 
Note  that  R3^  and  Rab  are  independent  of  gab  until  we  relate  gab  and 


R~bcd  “  RC  b  [cd] 


Rc 


Rc 


[bed] 


=  0 

=  0 


bed 


. 
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r using  (3.8).  R  measures  local  gravitational  effects  through  the 
geodesic  deviation  equation.  For  a  geodesic  vectorfield,  ka,  in  some 
locality,  representing  a  group  of  freely  falling  test  particles,  the 
distance  between  nearby  particles  is  described  by  a  vectorfield  pa  which 
is  formed  so  that  the  Lie  derivative 


of  p  along  integral  curves  of 


,  a  . 

k  is  zero : 


f  b  _ 
*°(k)  P  =  °  • 


t  b  c  a  „a  T  b  T  c  d 

kkp,  =  R  ,  .  k  k  p  . 
*  ;bc  bed  r 


(3-18) 

Then,  the  relative  acceleration  of  two  nearby  test  particles  is: 

(3.19) 

This  is  the  geodesic  deviation  equation  in  general  relativity. 

In  the  investigation  of  mechanics,  we  regularly  use  the  concept  of 
the  "frame".  An  observer  has  a  space  coordinate  system  set  up  around 
himself  and  measures  events  with  respect  to  a  clock  carried  with  him. 

This  space  coordinate  system  and  this  time  constitute  the  observer's 
frame.  Those  frames  in  which  local  space  appears  simplest  (fictitious 
forces  are  eliminated)  are  called  the  freely  falling  frames .  Geometrically 
a  freely  falling  frame  is  postulated  to  follow  a  timelike  geodesic  (in 
the  idealization  that  the  observer's  laboratory  is  of  point  size).  The 
are  length  between  two  points  on  the  geodesic  is  proportional  to  the 
elapsed  proper  time.  The  observer's  three  coordinate  basis  vectors,  which 
do  not  vary  as  viewed  by  the  observer,  are  three  mutually  orthogonal  space¬ 
like  vectors  all  orthogonal  to  the  tangent,  which  are  parallel  propagated 
along  the  geodesic.  At  any  spacetime  event  it  is  possible  to  set  up  a 
frame  in  which  the  (symmetric)  connection  is  zero  at  the  point.  Further, 
in  general  relativity  one  can  show  that  the  metric  has  the  flat  Minkowski 


form  there: 


^ab  ^  ab 


(3.20) 


. 


^7 


even  though  the  curvature  may  be  nonzero.  From  this  result  we  say  that 
the  physical  space  is  always  locally  flat . 

Symmetries  of  a  particular  solution  are  identified  by  Killing  vector- 
fields.  If  a  solution  has  some  continuous  symmetry  (eg.. ,  it  does  not  evolve 
in  time  (is  stationary)or  has  rotational  symmetries),  then  there  is  a 
vectorfield,  K  ,  such  that  an  observer  moving  along  its  integral  curves 
does  not  see  the  geometry  change.  Mathematically: 

^(k)  sab  =  0  •  (3.21) 

This  can  easily  be  rewritten  as  the  usual  form  of  Killing's  equation  in 
general  relativity: 


k(a;b)  =  0  ’  (3.22) 

Killing  vectorfields  generate  conserved  quantities  in  geodesic  motion. 

The  inner  product  of  the  tangent,  ta,  of  any  geodesic  with  any  Killing 

a 

vectorfield,  k  ,  that  exists  in  a  particular  solution,  is  constant  along 
the  geodesic: 


t°(g  ,  ka  t*3)  =  0  . 

ab  7  ;  c 


(3.23) 


a 

If  k  is  timelike  (a  stationary  solution) ,  the  inner  product  is  related 
to  the  total  energy  of  the  particle  in  the  field;  if  the  solution  is 
axially  symmetric  ,  the  angular  momentum  of  the  particle  is  related 
to  the  inner  product  of  the  tangent  and  a  killing  vector. 

Parallel  Transport  and  a  Real  Metric 

Let  us  begin  our  study  of  the  geometry  of  BMB  theory  with  considerat¬ 
ion  of  the  affine  connection  and  parallel  transport .  The  first  problem 
to  face  is  the  existence  of  several  connections,  only  one  of  which  (at 
most)  can  represent  the  parallel  transport  which  produces  physical 

effects.  There  is  La  ,  Pf  (the  metric  is  not  compatible  with  either  of 
' .  oc  be 

these)  and  there  is  the  Levi-Civita  connection,  now  labeled 


- 

' 
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compatible  with  Rewriting  eq.(3.9): 

^bc(lc)  ”2^  (S(ac),b  +  g(M)>c  -  S(bc),d)  (3.24) 

with  ?r^  defined  by  eq.(2.68).  This  connection  is  unrelated  to  L^c  or 
r"»ci  (34) 

1  be  ’  anc^  wou^-^-  need  to  be  given  serious  consideration  if  g^^x  were 
the  fundamental  geometric  tensor.  However,  since  it  seems  artificial 
to  separate  g(ab)  and  Sj-ahj  for  geometric'  purposes.,  nbc£lc)  is  not 
likely  the  "physical"  connection.  As  mentioned  in  discussing  Einstein's 
nonsymmetric  theory,  and  are  related  by  a  projective  transform¬ 

ation  given  by  eq.  (2.4.3^: 

rl  =  \ac +  !§*  . 

,  These  two  connections  therefore  have  the  same  geodesic  paths  (differently 
parametrized),  and  the  same  parallel  propagated  vector  directions,  but 
the  lengths  of  parallel  propagated  vectors  are  different .  Since  the 
simplest  set  of  field  equations  is  in  terms  of  T^c,  and  since  L^c 
does  not  appear  absolutely  necessary  to  derive  that  set,  seems  to  be 
the  most  likely  candidate  for  the  physical  connection.  These  arguments 
are  not  all  that  strong,  so  the  final  judgement  awaits  work  on  equations 
of  motion^^  ,  light  propagation  (see  chapter  5) »  or  other  physical 
properties . 

It  is  difficult  to  give  the  metric  tensor  full  geometric  interpret¬ 
ation,  primarily  because  it  is  not  compatible  with  the  connection: 


®ab ;  c  ^ad  ^[cbj  * 


(3.25) 


We  begin  with  the  use  of  a  real  gab  (and  will  come  to  hermitian  symmetry 
later)  and  immediately  must  decide  whether  gab  or  g^ab^  is  the  fundamental 
geometric  tensor.  With  S(ab)  i_b  is  easy  to  define  the  general  inner 
product,  but  such  separation  of  S^ab^  and  ^j-abj  as  contrary  to  unification, 
and  we  should  really  try  to  work  with  gab.  Having  made  this  choice  we 
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find  ambiguity  in  the  inner  product:  the  expected  expression  is  not 


symmetric  in  its  arguments: 

ka  ia  ^  Sat  la  .  ,  (3.26) 

For  two  arbitrary  vectors,  there  is  no  way  to  single  out  one  of  the  two 

possibilities.  The  ambiguity  only  disappears  for  magnitudes: 
a  T  b 


Sab  k  k  g(ab)  k  k  •  (3-2?) 

It  is  not  possible  for  the  general  inner  product  to  be  preserved  in 

parallel  transport  (nor  would  it  be  if  we  defined  it  using  g^-^, 

unless  we  also  use(^r^c^c)>  ignoring  T^c) .  If  we  parallel  propagate 

vectors  la,  ma  along  a  vector  ka,  the  inner  product  changes  according  to: 


T  c  /  ,  a  ck  _  nd  a  bT  c 

k  (sab  1  m  >,e  =  ^ad'cct]1  m  k  1 


(3.28) 

In  general  this  is  not  zero,  but  it  will  be  zero  in  the  special  case  of 

ci  £L 

m  being  k  .  This  means  that  the  length  of  the  tangent  to  a  geodesic 


is  preserved: 

;c 


kC  ka  =  0 


k°(g  ,  k3^^)  =  0 

Vtoab  ' ;  c 


(3.29) 

(3-30) 


Si 

and  one  of  the  two  possible  inner  products  of  a  vector,  1  ,  parallel 

Si 

transported  along  a  geodesic  with  the  tangent,  k  ,  is  preserved: 

k°(g,v  lakt).n  =  0 


5ab  "  '  ;  c 

k°(g  ,  kalb)  4  0  . 
v  ab  ;  c 

These  two  conserved  quantities  can  be  used  to  simplify  the  calculation  of 


(3.31) 

(3.32) 


particular  geodesics  and  associated  vectors  in  a  solution.  Note  that  the 
magnitude  of  a  parallel  propagated  non-tangent  is  not  usually  preserved. 
This  leads  to  an  unexpected  property :  under  appropriate  circumstances  a 
spacelike  vector*  parallel  propagated  becomes  null  (at  a  single  point), 
and  then  timelike(56t  To  set  up  a  frame,  therefore,  we  cannot  simply 
parallel  propagate  initially  chosen  spacelike  vectors.  Unfortunately  there 


-  .  * 


50 


is  no  alternate  to  parallel  transport  immediately  obvious,  so  we  do  not 
presently  know  how  to  construct  a  frame's  basis  vectors. 

There  is  another  problem  in  trying  to  construct  frames:  we  do  not 
seem  to  have  local  flatness.  Since  we  do  not  expect  to  be  able  to  eliminate 
the  electromagnetic  field  by  changing  coordinates  (if  gj.  ^  is  nonzero 
in  one  set  of  coordinates,  it  is  nonzero  in  all  of  them) , we  can  only  ask 
for  Mirikowskian  g^^^  at  a  point: 

Sab  =  ?  ab  +  S[ab3  *  (3-33) 

Then,  however,  eq.(2.2)  is  not  solved  by  zero  Y^c  .  In  fact,  is 

a  tensor,  so  it  cannot  be  reduced  to  zero  by  coordinate  transformations. 


If  we  look  at  coordinate  transformations  on  Y, 


a 


r(bc)  =  0 


a 


^(bc,d)  ° 


probably  the  strongest 

statement  that  can  be  made  is  that  by  appropriate  choice  of  coordinates 

one  can  make  the  symmetric  part  of  the  connection  and  all  totally 

(57) 

symmetrized  derivatives  vanish  at  a  single  point v  : 

a  ~  (3.34) 

(3-35) 

P/a  ,  x  =  0  .  (3.36) 

(be , d . . • e ) 

In  eq. (2.2) ,  however,  these  do  not  imply  that  has  constant  co¬ 

efficients.  In  addition  to  needing  a  method  to  construct  well  behaved 
frames  we  must  decide  whether  it  is  physically  reasonable  that  nonzero 
electromagnetic  field  means  space  is  not  locally  flat ,  even  as  seen  by  a 
neutral  test  particle.  Certainly  space  is  locally  flat  for  a  neutral  test 
particle  in  EMT. 

We  have  already  called  the  autoparallel  curves  geodesics  in  using  eq. 
(3*29)*  It  is  more  natural  to  call  the  generalized  "straight  line  the 
curve  whose  tangent  is  always  parallel  to  itself,  and  in  this  case 
abandon  extremization  of  arc  length.  Extremals  of  arc  length  still  exist, 


•i 
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solutions  of: 

ds  =  0  (3.37) 

ds2  =  gab  dxa  dxb  .  (3.38) 

It  is  easy  to  show  that  these  define  the  Levi-Civita  symmetric  connection 


on  the  manifold,  with  respect  to  which  these  curves  are  geodesics.  Thus 
the  Levi-Civita  connection  exists  along  with  rf  .  We  have  already  argued 
preference  for  P^c>  so  to  use  extremals  of  arc  length  as  free  fall  paths 
would  he  contradictory. 

The  curvatures  R  an<^  ^  are  still  given  by  eq.  (3.12)-(3.24) , 

but  R^,  symmetric  in  EMT,  now  has  both  symmetric  and  antisymmetric  parts. 

Ra,  ,  is  still  antisymmetric  in  c  and  d  but  the  identities  (3.16) , (3.17) 
bed 


become 


(58). 


R' 


a 


=  *  n? .  +  2  r 


[bcdl  **  1  [dc]  tbe] 

C  ] 


a 

[cb ,  d] 


(3.39) 


Ra.  r .  ,  =.-2  Ra.  ,r.  rf„-,  .  (3.40) 

b[cd;e]  bf[d  ecj 

The  geodesic  deviation  equation  is  modified  by  the  torsion.  Setting  it 

ci  ct 

up  in  the  same  way  as  before,  for  vectorfields  k  ,  p  : 

=  0 


ke  kb 


;c 

4k) pa  = 0 

*b  P%c  =  4cd  kb  k°  Pd+ 


b  d  c 
1  k  p 


ft  a 
1  0>cj 


,b  ,d  o 
k  k  p  . 


(3-41) 

(3.42) 

(3.43) 


id 


We  can  rewrite  eq . ( 3 * 43 )  several  ways,  none  of  "them  particularly  simple. 
On  the  basis  of  the  geodesic  deviation  equation  in  general  relativity  it 


can  be  shown  that  Rabcd  is  completely  determinable  by  measurement  of  the 
relative  acceleration  of  nearby  test  particles.  With  the  added  complicat¬ 
ion  in  its  replacement,  eq.(3.43).  it  W  be  difficult  to  separate  the 


•  .  ] 
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curvature  from  the  torsion  on  the  "basis  of  physical  measurements. 
Hermitian  Metric 

The  hermitian  formulation  involves  the  assumption  that  the  metric  is 
hermitian  symmetric: 


5ab  ®ba 


(3M) 

(the  bar  indicates  complex  conjugation)  .  This  means  as  real, 

imaginary.  The  easiest  way  to  work  with  this  is  to  simply  take  the 
theory's  new  constant  to  be  imaginary,  ie,  ik,  instead  of  real.  This 
formulation  gives  the  metric-  a  symmetry  again,  and  results  in  the 
unusual  "regular"  solution  (eq.2.18)  already  discussed,  but  the  complex- 
ification  that  necessarily  results  is  difficult  to  interpret . 


From  field  equation  (2.2),  relating  g  ^  and  P^c»  for  a  diagonal 

component  of  the  metric  g  (no  sum),  which  is  real: 

Si  Si 

-  rd  +  rd 

®aa,c  ~  ^ad  ca  ®da  ac 

—  nd  +  —  pd 

®ad  ca  °ad  ac 


(3.45) 


The  two  terms  on  the  right  hand  side  must  be  complex  conjugates.  This 
implies: 

pa  _  pa  (3*46) 

1  be  ~  1  cb  ’ 

the  connection  is  hermitian  in  its  lower  indices.  It  is  interesting  to 
note  that  the  other  connection  L^c  is  not  so  nicely  hermitian  symmetric. 

Si 

Using  the  projective  transformation  (2.43),  we  find  that  is  Pure 

imaginary,  but  that  generally  has  'both  real  and  imaginary  parts. 

Writing  out  the  geodesic  equation: 


k°  ka  =  k°  ka  +  Pa 

; c  ,C  be 

=  0 


(3.47) 


we  see  that  only  ri)  is  involved,  so  this  is  a  real  differential  equat¬ 
ion.  With  real  initial  data  it  will  generate  real  solutions  (geodesics  in 
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real  space) .  The  general  parallel  transport  equation: 

,  cT  a  ,  cT  a  ,  Pa  T  b  ,  c 

tk  =  t  k  +  p,  k  t 

; c  ,  c  "  he 

(3.48) 

=  0 

involves  the  full  and  thus  is  a  differential  equation  with  complex 

cl 

coefficients,  so  for  real  initial  data  for  k  we  must  expect  the  solution 

to  he  complex,  even  if  transport  is  along  a  real  t 

(39) 

Moffat  has  suggested  using  complex  conjugation  in  the  inner 
product  in  this  formulation  because  it  is  then  preserved  in  parallel 
transport  (under  certain  circumstances),  eg: 

®ab  k*  ?  '  <3'49) 

With  this  inner  product  the  magnitude  of  a  vector  is  always  real,  hut  the 

inner  product  of  two  vectors  is  complex,  and  ambiguous :  the  complex 

conjugate: 

(sab  k*  ^  =  sab  ta  kb  (3.50) 

is  just  as  likely  a  choice  (we  will  see  that  the  other  possibilities, 

g  ta  k°  and  its  conjugate,  do  not  mesh  with  parallel  propagation  as 
SLD 

nicely).  We  can  generalize  the  concept:  to  define  the  trace  of  a 
tensor^^  ^a. . .i  . . .  j .  . .  on  some  Upper  indices  i,  j,  (  there  is 

U  •  •  • 

similar  argument  for  two  lower)  we  must  resolve  the  tensor  into  vectors. 
There  exist  ,  .  . .  ,1^  ,  .  .  . , J .  .  ,B^,  .  .  .  such  that : 

¥AW-I(i)-J(n)-B(S)-  •  (3'51) 

Then  (being  careful  not  to  use  a  misleading  shorthand  notation  for  the 
trace) : 

a .  .  .  i  •  •  •  j .  . .  v-  «  a  Ti  7  j  p  ^  ^0 

gijQ  b...  =  sij  ¥A(n)-**  J(n)  J (n)  B  b  ^ 

Let  ua  parallel  propagate  ka  and  ta  (not  ta)  along  ma; 

n>°  ka  =0  (3-33) 

;c 

c  .  a  _ 

m  t  -  0 

;c 


(3-54) 
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a  -b 


a  rb 


;  c 

+  m°  ). 


ab;  c 


3ab  ;  c 


T  a/  C  rb  C— b  rid.  C  — b  _d 

=  k  (m  t  g  -  m  t  g  „  fr  -  m  tD  g  V 

ab,c  6ad  1  be  fedb  1  ac 

+  mC  g  t*3  +  mC  g  td  ) 

ab  ,c  6ab  1  dc' 

Taking  the  complex  conjugate  of  (3.54): 


-c  rt>  _  -c  rd.  nD 

mt  -  -  m  t  I  ,  . 

,  c  cd 

Q 

If  m  is  real,  then: 

C/  a  -bv  _  a,  c  rb  c  rb  nd 

m  (sab  k  t  );0  -  k  (m  t  gabjC  -  m  t  g(ft  Ta0 

c  -d  pb  \ 

'  m  Sab  *  "cd> 

=  k  t  m  ra„  -  b„,  ) 


-C  rd  nb 


3ab,c  6db  ac  6ad  1  cb 


=  0 


(3-55) 


(3.56) 


(3.57) 


by  virtue  of  the  field  equation  (2.2).  Thus  the  inner  product  is 
preserved  in  parallel  transport  along  real  vectors. 

The  major  problem  of  the  hermitian  theory  now  confronts  us.  Just 
what  does  the  increased  information  resulting  from  complexification 
represent?  What  is  the  significance  of  a  complex  valued  inner  product 
(complex  valued  angle?)?  Does  spacetime  now  have  four  complex  dimensions 
(equivalent  to  eight  real  dimensions),  so  that  complex  vectors  merely 
point  in  complex  directions?  If  so,  what  is  the  physical  significance 
of  the  increased  dimensionality,  only  detectable  when  there  is  an 
electromagnetic  field?  Note  that  parallel  transport  does  not  preserve 
inner  products  in  all  directions:  the  metric  is  once  again  not  compatible 
with  the  connection.  Is  it  possible  that  spacetime  still  has  just  four 
(real)  dimensions,  and  the  complex  part  of  a  vector  does  not  point 
somewhere",  but  represents  something  else?  (What?)  Given  a  complex 
vector,  how  do  we  parallel  transport  along  it,  or  find  the  geodesic 
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tangent  to  it?  If  we  simply  take  its  real  part  (itself  a  vector  in  this 
case) ,  we  find  the  length  of  the  real  part  is  not  the  same  as  that  of  the 
whole.  A  timelike  complex  vector  can  have  null  or  spacelike  real  part. 

Thus  the  (preserved  i  inner  product  loses  geometric  meaning.  Since  we  have 
no  resolutions  to  these  problems  at  the  present  time,  we  will  return  to 
working  with  a  real  metric . 

Killing  Vectorfields 

When  we  dicussed  Killing  vectorfields  in  general  relativity,  we 
identified  them  as  vectorfields  (eg  ka)  along  which  the  Lie  derivative 
of  the  metric  vanishes  (eq.(3.2l)): 

^(k)  S(ab)  =  0  ’ 

g(ai3)  representing  all  the  gravitational  information.  With  an  electro¬ 
magnetic  field  present  we  would  add: 

^(k)  Fab  =  0  ‘  (3-58) 

In  nonsymmetric  field  theories,  with  g  ^  representing  all  the  gravitational 

and  electromagnetic  information,  and  with  the  Lie  dervative  existing 

independent  of  metric  and  connection,  the  obvious  choice  is: 

(3.59) 

(53). 


■^(k)  sab  "  0  ' 


This  can  be  written  as 

~  T  C  T  C  .  T  c 

0  =  g  ,  k  +  g  k  ,  +  g  ,k 

6ab;C  •  b  cb  •; 


5ac  ;b  cb  ;a 


(3.60) 


=  e;  k+g  k  ,  +  g  k 

gab,c  &ac  ,b  6cb  ,a 


Using  the  field  equation  (2.2)  to  eliminate  derivatives  of  gab,  we  can 
rewrite  these  in  two  alternate  forms: 

0  =  gac  k°;b  +  Scb  k  ,a  +  ^ob  lad] 

c 


'c  ,ka 


(3.61) 


_  „  k°  +  K  ku  +  g  r  L  k  d  +  g  .  r° .  kd  . 

“  gac  ,b  Scb  ,a  6ac'  db  &cb  1  ad 


In  analogy  with  general  relativity  consider  the  value  of  the  inner 
product  of  ka  with  ta,  the  tangent  to  a  geodesic,  along  the  geodesic: 


: 
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(3.62) 


(3.63) 


Neither  choice  of  inner  product  is  preserved  nor  does  there  appear  to  be 
any  combination  of  them  which  is  preserved.  Defining  the  inner  product 


(ab)  ^oes  n0^  improve  the  situation  (unless  we  also  use  ^c) ) 


with  g 


Killing  vector  fields  do , not  appear  to  generate  constants  of  motion  in 
any  simple  manner. 

The  Viability  of  the  Theory 

In  this  chapter  we  have  seen  a  large  number  of  problems  presently 
associated  with  BMB  theory,  most  of  them  under  the  topic  of  finding 
physical  meaning  for  the  geometric  structure,  that  is,  translating  the 
mathematical  theory  into  the  physics  we  are  really  interested  in  (beyond 
the  minimal  physics.  Since  timelike,  null  and  spacelike  geodesics  exist, 
we  can  identify  paths  of  massive  and  massless  freely-falling  particles) . 
Until  most  of  these  problems  are  resolved  BMB  theory  will  be  at  a 
disadvantage  when  compared  to  EMT  because  we  do  not  know  how  the  math¬ 
ematics  translates  into  physical  prediction.  We  will  not  just  abandon 
the  theory,  because  there  are  some  areas  of  study  immediately  obvious 
which  may  provide  solutions  to  some  of  these  problems .  One  of  them  is 
geometric  optics:  by  investigating:  light  and  gravitational  radiation 
propagation  we  should  be  able  to  see  which  connection  really  is  used  by 
null  fields.  This  will  be  discussed  in  Chapter  5. 
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Chapter  IV 

Perturhation  Expansions 
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Perturbations  in  EMT 


Perturbation  theory  is  a  well  known  approximation  method  in  physics, 
including  general  relativity.  On  some  initial  configuration  (assumed 
known) ,  we  consider  a  change  in  the  field  quantities  which  is  small  enough 
to  produce  only  small  changes  in  all  other  relevant  quantities.  Then 
the  field  equations  are  expanded  in  the  small  amplitude  of  the  change,  X  , 
and  each  order  of  "smallness"  is  solved  individually.  With  a  linear  theory 
we  can  formally  write  out  a  power  series  in  A  and  declare  that  there  is 
some  finite  radius  of  convergence  in  A.  With  a  nonlinear  theory  (both  EMT 
and  BMB  are  nonlinear),  convergence  is  harder  to  assess.  Isaacson^ ^ 
indicates  that  possibly,  only  the  first  (linear)  correction  is  to  be 
trusted. 

We  will  write  out  expressions  for  EMT  perturbations  here  (this  is  all 
fairly  well  known  material)  so  that  we  can  compare  them  with  their  BMB 
analogues.  There  are  three  cases  we  will  be  concerned  with,  in  order  of 
increasing  complexity:  empty  Minkowskd  space  background,  gravitational 
but  no  electromagnetic  field  background,  and 'general  background  fields. 

Formally,  we  assume  all  relevant  quantities  can  be  expanded  into  (at 
least)  a  polynomial  inX  with  possible  other  terms  in  higher  orders: 

(4.1) 


gab  =  +  +  ^  S'n  + 


3ab 


3ab 


3ab 


r a  =  a  +  a  r  ?  +  p  +  •  •  • 

1  be  I  be  Al  be  **  1  be 

Fab  =  ••• 

This  implies,  for  example: 

Rab  =  \b  +  XRab  +  ^  R  ab  +  •  •  •  ' 

'A is  an  artificial,  dimensionless  "small  parameter  : 

"X«l 


(4.2) 

(4.3) 

(4.4) 

(4.5) 


which  is  in  fact  eventually  ignored,  but  it  is  convenient  when  handling 
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several  independent,  types  of  smallness  (each  can  be  labelled  differently) , 
as  will  be  necessary  later.  We  take  all  quantities  in  the  equations  other 
than  ^  as  being  of  order  unity.  The  order  zero  (in  A)  quantities 

are  assumed  to  be  those  of  some  exact  solution  in  EMT  (later,  in  BMB 
theory) ) 

With  empty  Minkowski  background  (first  order  will  be  lineaurized  theory) : 


=  n. 


ab  '(■  ab 
=  0  . 


ab 

Necessarily, 

o 

F  =  0 
ab 

and,  for  simplicity,  we  choose  the  obvious  coordinate  condition: 

rbc  =  0  • 

The  general  EMT  field  equations  are: 

R  ,  =  8ttT 


ab 


ab 


F  =0 

Cab ,  c] 

Fab  ,  =  0 

;b 


(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 


T  .  =  F  F.+fg.F.F 
ab  a  cb  ^  ab  cd 


cd 


(4.13) 

Applying  our  expansions  in  this  case,  we  find  (after  noting  that  the 

2 

zeroth  order  is  trivial)  that  T  is  of  order  X ,  being  quadratic  in  F  , 

cl  D  clD 

and  that  to  first  order  covariant  differentiation  is  coordinate 
differentiation.  Then  the  first  order  gravitational  and  electromagnetic 


equations  decouple: 

r  be  =  *  ?  (gbd,c  +  Sdc,b 

-  r c  -  = 0 


Sbc,d^ 


ab 

Cab ,  cQ 
'ab 

,b 


ab ,  c 
=  0 

=  0  , 


ac ,  b  i 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


. 


6o 


and.  the  first  order  gravitational  and  electromagnetic  perturbations 
are  independent  of  each  other.  If  we  go  to  second  order  (in  spite  of 

?1  »  ft 

Issacson's  warning)  we  see  that  is  modified  by  F  ^  through  T  and 

ii  i  i 

F  ,  is  modified  by  g^  through  P  ac  appearing  in  covariant  differentia- 


ab 
tion. 


In  the  case  of  a  purely  gravitational  background  we  have  nontrivial 

O  0  3,  O  O 

Sab’  P  and  R^  but  retain  zero  F^.  Now,  order  zero  is  not  trivial 


but  can  be  anticipated: 
o 

R  ,  =  0  . 
ab 

The  first  order  equations  (T  ^  is  again  second  order)  are: 


i — i  si  _  i  ^ad 

be 


r— i  a-  _  1  (  O 

I  -Krt  2  S  (  S 


+  °  0  \ 

bd,c  ^dc,b  ^bc.d' 


i  oad  /  ’ 

+  2  S  ( -§bd,c  +  ^dc,b  ®bc,d) 


ab 
o  c 


=  r  c  -  n  c 
1  ab,c  1  ac,b 


-  I*  ad  r  Ob  -  I* 


d  A  c 

cb  ad 


rac  r 


’  d 


ab 


o  d  p  c 
'  ab  dc 


(4.18) 


(4.19) 


(4.20) 


oab 


=  0 


Fab  =0 

,  o 

F  ,  =0 

Qab ,  c] 


g  is  the  inverse  of  g^ 


'  ab 

g  is  defined  by  the  expression: 


ab  o  ab 
S  =  S 


^  'ab 
+Ag  + 


and  can  be  written: 
'ab 


o  ac  o  bd  1 

-  e  s  scd  • 


(4.21) 

(4.22) 

(4.23) 


(4.24) 

In  eq.  (4.21),  the  notation  means  covariant  derivative  with  respect 

o 

to  a  ,  as  the  divergence  equation  is  only  taken  to  first  order  (  the 
I  be 

indices  have  been  raised  with  g  ab)  .  The  gravitational  .and  electromagnetic 
disturbances  here  are  again  independent  to  first  order,  and  in  second 

V  n 

order  R  ^  is  modified  by  F^  through  T^ . 


■ 


. 
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Finally,  in  the  case  of  a  fully  general  background  we  have  the  stand¬ 


ard  EMT  equations  for  order  zero: 
o  o 

R  ,  =  87 r  T  , 
ab  ab 


=  0 


o 

F 


[ab ,  c] 


0 


and  T  ^  has  a  first  order  part: 


ab 


'  cd  0  0  0  c  » 

=  g  F  F  ,  +  F  F  ,  +  F°  F 

ac  db  a  cb  b  ac 


°c  ’ 


(4.2  5) 

(4.26) 

(4.27) 

(4.28) 


+  *  ®ab  <**  L  +  K  Fed>  +  F 


2cd 


cd 


where  indices  of  have  been  raised  with  gab.  Here  gab  still  has  the 


(4.29) 


(4.30) 


form  of  eq.  (4.24),  and: 

Fab  =  gac  gbd 

=  Fab  +  Xf3-13  +.  . . 

'ab  oac  obd  '  ,  oac  'bd  °  'ac  obd  °  . 

F  =  s  8  Fcd  +  ®  s  Fcd  +  s  8  Fcd- 

V 

Q 

Rab  has  the  same  form  in  r'bc  as  was  seen  in  The  first  order 

field  equations  are  eq.  (4.19)  and: 

Rab  =  8lrTab  (^l) 

•ab  ^  +  jlac  o  +  jjcb  ^  a^  +  °ac  p  ^  +  °cb  p  ^  =  Q  (i,.32) 


»*> 


Fr  ,  _  =  0 

Qab,c] 


cb 


(4.33) 


lab 


The  first  three  terms  of  eq(4.32)  are  just  F  .,  .  The  detailed  form  of 

6D 

these  equations  is  rather  complicated  and  difficult  to  solve  at  present . 

•  • 

Note,  however,  that  a  g  ^  usually  generates  an  F^,  and  vice  versa. 
Perturbations  In  BMB 


-  Minkowski  Background 

For  Minkowski  background  in  BMB  we  take ,  in  analogy  with  EMT : 


g 


ab 


I  £  ,f 

=  ;Zab+lsab+  Xgab+-" 


(4.34) 
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where  ^  Is  fhe  flat  (symmetric)  Minkowski  metric.  Also: 

Sab*  ?ab  +  Xgab  +  x2gab+  ...  , 

and  by  solving: 

ah  _  g  a 


S  g 


cb 


we  find: 


'ab  ^,ac  -.db  ' 

g  =  -  n  n  sd0 

"ab  ac  db  "  ac  db  ef  '  ' 

g  =  -  n  n  sd0  +  n  n  7  ge0  sdf  • 


Now: 


. .. .  Si  -  ,  3/ 

rb0“  xrbc  + 


-2  "  a 

^  r  b0+  ••• 


and.  (4.34)  in  field  eq .  (2.2)  relating  g  and  P  imply: 

n  a  i  ad  '  '  , 

I  (be)  2  4  ^(dc) ,b  ^(bd),c  ®(bc),d' 


(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 


(4.40) 


^(bc)  2^  ^®(dc),b  +  ^(bd),c  ^(bc) ,d  2 


|  ■  0 

[ec]  '[db] 


2  g  rv„T  r 


[be]  1  Cod]  2  ®(de)  1  (be) 


n  a  _  l  ~ad  /' 

1  [be]  -  2?Z  (s 


[dc],b  +  ^[bd],c  +  ^Q:c],d^ 

rrbci=  *2  (  srj„i  v  +  g 


(4.41) 

(4.42) 


a  _  _i  ad 

[be]  2t  ^  s[dc],b  T  6[bd3,c  '  6 [be]  , d 


If 

+  Si 


2  S[ec]  ^  (db) 


o'  r»  e  o'  n e  ^ 

g'[bel  '  (cd)  g(de)  [be] 


(4.43) 


R  ^  in  terms  of  ff  is: 

be 


ab 


R(ab)  ^  (ab),c  2  ^^(ac),b 


+  r(be),a) 


(4.44) 


4  d 


E(ab)=  r(L).o  -^(r(:e).b  +  r(bBc),a)-  r(ad)  r(,b) 


A  c  Aed  ,  n  a  n u 

"  1  Cadi  1  [cb]  1  (cd)  1  (ab) 


4  d 


R[ab]  n 


[ab]  ,  c 


I! 

R 


4  d 


Ac  Ad 


[ab] 

-  hi,  r 


nC  pC  pd  popu. 

1  [ab],e  '  1  (ad)  1  [cb]  '  1  [ad]  1  (eb) 


(4.45) 


(4.46) 


(4.4?) 


(cd)  [ab] 


These  last  expressions  are  written  using  the  tracelessness  of  the  torsion 


. 


' 
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(which  here  becomes  tracelessness  in  each  order) .  I  ,  in  berms  of 

ab 


gab  13 5 


1/  ,  %  =  0 
(ab) 


(ab) 


[ab] 


lab] 


-8tt  /  cd  ' 

72  ^  8 

k 

- 8rr  ' 

k2  gCat: 

-8tt  " 

2  g[ab]  ’ 


[ad]  g  [cb] 


+  kri 


ab  g[cd] 


'  Ted]  x 

g  ) 


O  g 
k 


(4.48) 

(4.49) 

(4.50) 

(4.51) 


From  these  the  field  equations  in  first  order  (after  solving  for 


r*0)  are: 


?b°  *M.c  =  ° 


1  n2  l 

2  □  g 


[ab]] 


8tt  ' 

k2  g[ab] 


=  V  -  V 
a  i  b  b  ,a 


*~^g(ab)  +  ^  ^  g(cd),ab  +  g(db),ac  +  g(ca),db^ 

Note  that  in  first  order,  the  gravitational  and  electromagnetic 


(4.52) 

(4.53) 

(4.54) 


perturbations  (g^-^  and  S^-^)  are  independent.  Eq.(4.54)  is  one  that 
can  be  found  in  EMT,  meaning  first  order  gravitational  perturbations  are 
the  same  as  EMT's.  The  second  order  general  case  equations  are  complicated 

f 

but  we  can  look  at  special  cases  (see  Appendix  ).  If  we  set  gj-^-j  bo 
zero,  the  gravitational  and  electromagnetic  equations  become  independent 
down  to  second  order,  and  the  second  order  gravitational  equations  are 
equivalent  to  the  analogous  equations  in  EMT.  The  gravitational  field 
does  not  generate  an  electromagnetic  field  unless  there  is  an  electromag- 

f 

netic  background  to  perturb  (the  same  is  true  in  EMT) .  If  we  set  g^-^ 

to  zero,  then  the  second  order  symmetric  field  is  governed  by  the 

2 

equations  of  EMT  plus  terms  of  order  k  : 


R  ,  (  ,  \  =  87rT  ,  /  +  \  +  0(k  ) 

ab  (  emt )  ab  emt ) 


(^•55) 

The  most  interesting  solutions  in  this  approximation  are  plane  waves, 
representing  radiation.  They  will  be  discussed  in  chapter  5* 


. 


. 
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-  Vacuum  Gravitational  Background 


For  perturbations  on  a  purely  gravitational  background: 
o 


toab 

rt  = 


(ab)  +^ab  +  *  * 


be  +  * '  * 


(4.56) 

(4.57) 

(4.58) 

Because  of  increasing  complication,  we  will  not  look  at  the  second  order 


a 

(be) 


Rab  R(ab)  +^ab  +,,< 


in  anything  but  the  Minkowski  background  case.  BMB  reduces  to  EMT  when 
there  is  no  g  >  so  we  can  correctly  predict  that  the  zeroth  order 
solution  is  a  vacuum  general  relativity  geometry.  Inspection  of  the 
first  order  equations  shows  that,  in  analogy  with  EMT,  the  gravitational 

1  » 

(g(ab))  and  electromagnetic  (gj-a^]  )  fields  are  completely  independent  in 
first  order,  so  we  can  consider  them  separately .  The  equations  in  g/  ^ 
are  equivalent  to  those  found  in  EMT  (in  fact  general  relativity 

I; 

without  F  ^),  and  we  need  not  write  them  out  again.  In  g^-j  we  find: 


pa  _  1  g(ad)  / 

I  2  S  v&i 


D>c]  "  2  0  VB[bdJ,c  +  g[dc],t 

"  2  gCbe]  ^(od)  _2gCeol  ^(db)'  ' 
The  equation  demanding  zero  trace: 

I'm  = 0 


+  g 


Cbc3 » d 


(4.59) 


is  then: 


o  (cd)  /  '  '  p  e  _  1  p  e  Y=  0 

s  ^  g[ad],c  "  S[ae]  '  (dc)  ~  S  [edj  ‘  (ac)  ^ 


(4.60) 


(4.61) 


The  part  in  brackets  is  just  g£acg.c  S0  as  a  sPecial  statement  of 


the  divergence  equation  on  g 


index  instead  of  using  g 


tab] 


Cab] 


mentioned  before,  this  time  raising  an 


For  I  ,  we  find  simply: 
ab 


I,  =  ° 

(ab) 


'Cab] 


8t r 

k2 


g 


Cab]  * 


(4.62) 

(4.63) 


The  equation  in  $  [ap]  can  wrlbben: 


' 


■ 
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d 

£caU 


+  rv-.v  i  v\,  -  %  P-.) 


[^(cd)  rt 


;*■  P£ab] ,  o  T(ad)  6  £cb]+  5(bd)  P [ 

(4.64) 

Cab]  "  k2  6tab]',e]  =  0  ’ 

but  in  terms  of  g  ^  alone  it  is  complicated.  These  equations  are  used 

later  in  this  chapter  in  finding  BMB  perturbations  of  the  Schwarzschild 

geometry,  and  in  chapter  5* 

-  General  Background 

In  the  general  case  the  background  g  ^  and  p  have  symmetric  and 

V 

antisymmetric  parts.  It  is  no  longer  easy  to  solve  for  I  ^  in  the 

»  a  ’a 

straightforward  manner  used  earlier,  as  both  P  (bo)  and  pc]  appear 


in  each  of  the  equations  on  g 


(ab) 


and  Sj-g^-j  (from  eq.(2.2)): 


n_  jQ.  d  o  i— ,  d  £d 

g(ab),c  g(ad)  I  (be)  g(ad)  '  (be)  g(db)  '  (ac) 

o  '  d  '  .o  d  o  4d  '  od 

"  S(db)  '  (ac)  g[ad]  rCcbj  S[Iad5  1  M  gpbj  Ipac] 


(4.65) 


O  n  d 

gCdbI|  1  [acH 


'  '  0  d  o  '  d  '  o  d 

°  -  SCab],o  '  gTad]  Rbo)  "  6 [ad]  '  (be)  s[db]  l(a0) 

-  Su^  I21, 


o  1  d  Ad  O 

"SCdb]P(ao)  S/‘J'  S 


(ad)  '  []cb] 


(ad)  ^[cbl  g(db) 


A  d. 

rCacG 


(4.66) 


o  '  d 

g(db)  1  facj 

Of  course  we  have: 

'  c  _ 

Pc]  =  0  • 

R  ,  can  be  written: 
ab  ( 

R(ab)  =  P(ab),c  ’  r  ^  P(ac)  ,b  +  P(bc),i 


(4.6?) 


)  -  ft 


c 

(ad) 


r, 


(cb) 


o  d 


R 


Labi 


A  c  o  a  ,  o  c  p  ^ 

^(ad)  Rcb)  M [lad]  1  [ZcbH  1  lladj  I  [cb] 

+  r(ab)P(cd)  +  p(ab)  fyed) 

=  P fab] ,  c  6(ad)  peb]-  P(ad)  Pd]  P(cb) 


(4.68) 


|>b] 

-  Head]  l*(cb)+  pb]  P (cd)  +  TCab]ft(cd)  * 


(4.69) 


'  '  ••  ■■  Y  10 

j  - 
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^(ab)  anc^  ''’[ab]’  w^en  written  out  in  terms  of  our  expansion,  are  complicated 
expressions,  "because  most  of  1^  is  cubic  in  g^,  and  includes  the  upper 
indexed  metric: 

S  =  “  S  g  gdc  (4.70) 

(the  symmetric  and  skew  parts  of  this  both  have  several  terms) .  The 
field  equations  in  contracted  curvature  are  then  very  complicated.  We 
will  not  actually  deal  with  them  until  chapter  5»  when  the  geometric 
optics  approximation  allows  further  simplification. 

BMB  Perturbation  of  the  Schwarzschild  Geometry. 

In  EMT,  consider  perturbing  the  Schwarzschild  geometry  with  an 
electromagnetic  field  (to  first  order,  any  gravitational  perturbation 
is  independent,  so  we  can  ignore  it.  We  are  particularly  interested  in 
a  possible  BMB  magnetic  monopole  in  this  section) .  We  need  only  solve 
the  electromagnetic  equations  (4.21),  (4.22),  and  it  is  easy  to  show 
that  they  are  satisfied  by  either  a  radial  electric  field: 

E  =4  (4.71) 

r  2 
r 

and/or  a  radial  magnetic  field: 


B  =  4  (4.72) 

r  2 
r 

(we  are  not  ruling  out  EMT  magnetic  charge  a  priori  here).  Thus,  the 
electric  and  magnetic  Reissner-Nordstrom  solution  (eq.  (2.62),  (2.63)  ) 
is  found  in  approximation.  By  doing  analogous  work  in  BMB  theory  we 
may  find  the  known  electic  monopole,  and  some  clue  to  a  magnetic 
monopole,  in  approximation.  Recall  that  in  chapter  2,  we  said  that  no 
magnetic  monopole  solution  has  been  found  in  the  theory.  If  it  has  no 


such  manopoles ,  then  it  would  have  an  advantage  over  EMT,  since  EMT  does 
not  really  preclude  monopoles  yet  none  have  been  found  experimentally. 
The  methods  here  constitute  an  independent  test  for  monopoles.  Any 


. 
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approximate  forms  found  will  give  indication  of  forms  to  look  for  in 
further  exact  calculations. 

In  BMB  let  us  begin  with  a  strictly  radial  electric  field  perturbation: 

g[10]  =  &  E(r)  (4.73) 

in  the  usual  (exterior)  spherical  coordinates.  The  background, 
Schwarzschild' s  solution,  is  well  known: 

as2  =  (1  -  Is)  at2  -  ^  2ro)  -  r2(<»2  +  sin2eajif2)  (4.7*0 

x 

and  we  will  not  write  out  the  expressions  for  its  geometric  quantities. 
In  straightforward  solution,  eq.  (4.59)  implies: 


rm  - (1  -  r> 


2m  \  dE 
dr 


2m  > 


r-i  ^  _  f± 

1  [02]  ^r  2 


E 


=  rr 


3 


[03]  ’ 

Then  eq.  (4.60)  implies: 

+^=o, 

dr  r 

the  solution  of  which  is: 


(4.75) 

(4.76) 


(4.77) 


E  =  ^ 
r 


(4.78) 


(Q  a  constant  of  integration).  Eq.  (4.64)  is  then  satisfied  identically. 


For  the  case  of  a  radial  magnetic  field: 


>[23] 


=  B(r)  r  sin  9 


Solving  eq.  (4.59)  we  find: 

pi  _  _i  (i  _  2ni) 

I  j-23-,  "  2  U  r  ) 

'2  _  ±_  bB  •  q 

bmi  -  "2  ^  sln  y 


—  r2sin  9  -  2Br  sin  9 

dr 


[13] 

P  3 
1  [12] 


dB 


(4.79) 

(4.80) 

(4.81) 

(4.82) 


2  sin  9  dr 

Eq .  (4.60)  is  trivial,  as  the  torsion's  trace  is  identically  zero.  In 
eq.  (4.64),  let  us  condider  the  various  combinations  of  the  three  free 


indices  in: 
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»  w 

^Catf)  +  ICab]^  ,o  (4.83) 

"before  anti symmetrizing .  It  is  easy  to  verify  that  if  any  index  is 

zero,  the  expression  is  identically  zero.  Then  if  the  last  index  is 

three,  the  expression  is  zero,  as  the  term  to  "be  differentiated  is 

independent  of  This  leaves  two  independent  cases,  (a,b,c) 

equal  (3,1,2)  or  (2,3,1)  .  In  the  first  case: 

r-h.  2  o  2  '3 


^R[31]+  I[31^,2  ^  r [313.2  ^ (33)  Hl21J 

-  1^(32)  ^[31]  +  ^(23)  rpi]),2 


(  (i  H  sin  9) 


,2 


-  (-  sin  9  cos  9)  ( 


=* _  M)  ) 

2  sin  9  dr'  '  ,2 


(4.84) 


=  0 


and  in  the  second: 


(E[23]+  IC23^I1'“  (  ^*£23J 1 1  4(21)  4:23]  '  4(22)  4  £13] 

+  4(31)  432]  +  4(33)  4?2]+  4(4)  r(23Q'-ffW.l 


"  ~  (  4^  r2sin  9  -  2  Br  sin  9) 
r  dr  ' 

2 

-  i  (1-  — )  (^4  r2  sin  9  -  2  B  sin  9) 

r  dr 

-  (1  -  T}  d?  r  Sln  6  -  ^  Br  sln  ®J,i 

d2B 


(4.85) 


^  2 
/  x  s  fA  2nu  ,  dB  /  \  ,  -n  ( *  8irr  v 

2  (  -  "2")  (1'  T)  dr  (m-r)  +  B  (I"  —  ) 


,1 


sin  9  . 


_dr  k 

Since  there  is  only  one  nonzero  term  in  eq.(4.83),  antisymmetrization 

over  the  three  indices  has  no  effect,  and  eq.(4.64)  simply  means  eq.(4.85) 


is  zero: 


^  C  \)  (1-  f )  +  §  (m-r)  +  B  ( 1  —  )  =  K  . 


(4.86) 


dr  k 

K  is  an  integration  constant;  to  interpret  it  we  take  the  small  k  limit 


of  eq. (4.86) : 

87ir2 


B  =  K 


(4.87) 


« 
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B  =  - 


K  k* 
8irr2 


(4.88) 


K  =  - 


87rL 


(4.89) 


Assuming  a  power  series  in  r  that  satisfies  the  boundary  condition: 


lim  B  “  0  f 
r-^0® 


(4.90) 


that  is: 

B(r)  =  E 


OG 


_n 

n 


(4.91) 


n  -  1 


we  find: 


B(r)  =  \  +  + 


...  i 


where  for  the  higher  order  terms  (  n^5): 
2  r  2 


an  '  8tt 


(  -  f  +  f  -  2) 


a  „  +  (n-3)  m  a  , 
n-2  n-3 


(4.92) 


(4.93) 


The  series  does  not  terminate  to  form  a  polynomial . 

Several  comments  can  be  made  about  the  series  solution.  Since  it  is 
an  infinite  series,  it  is  not  among  the  family  of  polynomial  solutions 
tested  for  in  exact  calculation.  The  form  of  the  recursion  relation 
is  complicated,  so  that  the  series  is  not  recognizable  as  any  familiar 
function,  and  we  must  worry  about  convergence.  The  fact  that  there  are 
three  parameters  (k  ,  m,r)  whose  relative  sizes  are  important  and 
the  complicated  form  of  eq.(4.93)  combine  to  make  convergence  very 
difficult  to  assess.  We  can  note  two  limits.  When  k  is  zero,  we  get 
EMT,  and  just  the  expected  first  term.  Also,  all  but  the  first  term 
include  a  positive  power  of  m,  so  that  when  m  is  zero  we  also  find  simply 
the  first  term.  Only  in  these  two  limits  do  the  electric  and  magnetic 
monopole  have  the  same  form(in  this  approximation)  .  The  m  =  0  limit 
cannot  be  expected  to  be  a  physically  interesting  case.  Experimentally, 
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electric  charge  is  always  accompanied  by  mass,  yet  in  EMT's  Reissner- 
Nordstrom  solution,  m,  Q  and  L  are  completely  independent.  The  massless 
charged  Reissner-Nordstrom  geometry  is  a  naked  singularity  which  is 
repulsive  for  neutral  test  masses.  Only  in  the  hermitian  formulation 
of  BMB  might  this  case  be  physically  tolerable .  An  infinite  series  magnetic 
field  is  difficult  to  work  with  in  exact  calculation  for  this  case,  so 
the  general  form  found  here  may  be  difficult  to  actually  verify  or 
discredit.  On  the  other  hand,  the  simple  form  of  the  massless  charge 
result  suggests  that  one  can  test  for  it  in  exact  calculation. 


Chapter  V 

Radiation  In  The  Geometric 


Optics  Approximation 
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Geometric  Optics  in  EMT 

-  Introduction  :  Simple  Gases 

Elementary  applications  of  geometric  optics  in  general  relativity 
are  well  known,  and  the  more  exotic  cases  we  will  eventually  come  to  have 
“been  discussed  in  the  literature^'''  .  By  "geometric  optics 
approximation"  we  mean  the  assumption  of  a  periodic  disturbance  of 
high  frequency,  and  expansion  of  the  field  equations  into  power  series 
in  the  wavelength.  In  complicated  cases,  this  is  usually  coupled  with  a 
weak  field  assumption  to  obtain  equations  simple  enough  to  solve  easily. 

It  is  not  always  necessary  to  assume  high  frequency  to  find  a 

(64) 

reasonably  physical  radiative  field  ' .  The  linearized  EM1  equations 
discussed  in  chapter  4  (eq . (4.14)-(4.16) )  are  well  known  to  have  plane 
wave  solutions,  gravitational  and/or  electromagnetic,  which  travel  along 
null  straight  lines  in  Minkowski  space. 

We  can,  on  the  other  hand,  use  high  frequency  without  small  amplitude 
as  follows^)  .  We  assume  there  exists  a  geometry  which,  in  some  region, 
is  a  solution  of  the  vacuum  in  EMT  equations,  eq. (4.10)-(4.13) »  where 


F  ,  has  the  form: 
ab 


Fab  =  fab  8XP  (f)  +  co 


(5.D 


(where  "cc"  indicates  the  complex  conjugate  of  the  preceding  expression), 

that  is,  the  electromagnetic  field  is  strictly  high  frequency  in  the 

region.  We  assume  that  the  Riemann  tensor  is  bounded  in  the  region. 

f  itself  is  slowly  varying,  but  the  exponential  is  of  high  frequency, 
ab 

£  is  the  small  parameter  proportional  to  the  reciprocal  of  frequency  ) : 

£  «  1  (5.2) 

in  which  expansions  are  made,  and  s  is  the  "phase  factor"  whose  derivative: 

k  =  s 
a  ;a 

=  s 


(5-3) 


' 
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is  the  propagation  vector.  Note  that,  as  generated  in  eq.(5. 3),  the 
propagation  vector  is  lower  indexed,  and  is  technically  a  "one-form" 
rather  than  a  vector  (upper  indexed)  pointing  in  four  space  .  The 
index  must  be  raised  with  g  to  make  it  the  propagation  vector(this 
distinction  becomes  significant  in  nonsymmetric  theories) . 

The  rapidly  varying  field  here  is  in  the  standard  form  we  will  use 
throughout  our  study  of  geometric  optics.  It  may  be  thought  of  in  a 
Fourier  analysis  sense  as  a  "bump"  in  the  amplitude  versus  frequency 
function  with  central  frequency  ^  (but  rigorous  Fourier  analysis  is  diff¬ 
icult  to  do  in  curved  space) .  With  nonlinear  equations  there  tends  to  be 
interaction  between  different  frequencies,  and  this  is  likely  to  be  a 
significant  interaction  when  there  are  both  gravitational  and  electro¬ 
magnetic  waves  in  a  region  (if  they  are  generated  independently  they  are 
unlikely  to  have  identical  frequency) ,  yet  we  will  consistently  use 
notation  that  considers  only  a  single  high  frequency  (^)  band.  This  is 
because  we  are  using  the  geometric  optics  approximation  to  produce  cases 
where  the  field  equations  are  (approximately)  linear,  so  that  the  inter¬ 
action  between  different  frequencies  should  not  then  be  significant .  In 
the  cases  in  which  we  do  find  easy  to  solve  equations,  the  results  should 
be  valid.  In  cases  in  which  the  approximate  equations  are  complicated, 
we  can  conclude  nothing  beyond  the  fact  that  this  approximation  failed 
to  linearize  the  situation. 

Now,  let  us  consider  what  form  f  b  must  have.  The  "curl"  equation, 
eq . (4. 11 ) ,  implies  the  existence  of  a  vector  potential,  A^,  which  is  high 

frequency: 

.  /iss 

Ab =  ab exp  00 

=  j  £c°-b  exp  (if) +  °c> 


(5.4) 


(^"b  is  the  P°larization  vector)  so  that: 


Fab  lt0[a  V|  +  ^a.b]  ’  (5-5) 

Note  that  a^  is  first  order  in  E(smallness) ,  but  it  generates  a  unity 

order  term  in  .  Differentiation  of  a  high  frequency  term  of  the  nth 

order  m  times  generates  terms  as  large  as  m,  in  general.  Now,  raising 
3/"b 

indices  with  g  ,  we  find  that  the  divergence  equation,  eq.(4.12),  has 

1 

a  term  of  order  ^  : 

^jb  =  (exP(lf)  Oy)  tol"ak1:)-,kb  +  cc)  +  0(1)  (5.6) 

which  (if  nonzero)  will  be  the  dominant  term  in  the  high  frequency  limit 
(ie,  as  £  goes  to  zero),  and  which  therefore  must  be  zero.  This  implies: 

^  =  0  (5.7) 

(5-8) 

The  propagation  vector  is  null:  the  high  frequency  disturbance  propagates 
along  null  paths .  Also ,  the  polarization  vector  is  orthogonal  to  the 
propagation  vector. 

Throughout  this  chapter  we  will  be  primarily  concerned  only  with 
demonstrating  the  first  property  expected  of  "massless"  radiation  fields, 
that  is,  propagation  along  null  geodesics.  In  general  relativity  we  can 
use  a  small. (but  important  for  the  geometric  optics  approximation) 
theorem:  if  the  gradient  of  a  scalar  (eq.(5»3))  is  null,  then  it  is 
geodesic: 

kb  k  ,  =  0  .  (5-9) 

a;b 

clTd 

The  proof  is  easy,  but  it  relies  on  the  covariant  derivative  of  g  being 
zero  and  on  the  first  two  covariant  derivatives  of  a  scalar  commuting 
(which  is  true  with  a  symmetric  connection.  Neither  of  these  properties 
is  found  in  nonsymmetric  theories) .  Thus  we  see  that  in  an  exact  vacuum 
solution  (region)  with  high  frequency  electromagnetic  field,  this  radiative 


« 
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field,  travels  on  null  geodesics  in  the  high  frequency  limit,  and  we  are 
justified  in  identifying  null  cones  with  the  light  cones,  in  this  type 
of  situation. 

-Gravitational  Background 

There  is  a  slightly  different  problem  than  the  one  just  discussed 

which  is  important  in  its  own  right .  When  we  investigate  the  properties 

of  some  exact  EMT  (or  BMB<)  solution,  we  look  at  paths  of  test  particles 

and  light,  saying  that  we  expect  neutral  test  particles  to  follow  timelike 

geodesics,  and  light  to  follow  null  geodesics,  of  the  given  solution. 

These  expectations  can  only  be  trusted  when  supported  by  some  calculations 

of  how  small,  but  not  quite  trivial,  masses  and  small  electromagnetic 

disturbances  really  propagate  in  such  circumstances.  We  will  not  discuss 

the  massive  case  here,  since  it  is  part  of  the  theory  of  equations  of 

(55) 

motion,  which  we  have  not  detailed  ' ,  but  we  are  in  a  position  to 
discuss  light  propagation  (in  fact  the  framework  here  puts  electromagnetic 
and  gravitational  disturbances  on  an  equal  footing,  so  we  can  discuss 
gravitational  radiation,  too). 

Consider  a  physical  light  signal  moving  through  some  background, 
as  a  perturbation.  If  its  amplitude  is  of  order  unity  (as  defined  by 
the  background) ,  then  its  energy  density  will  rival  that  of  the  back¬ 
ground,  causing  significant  changes  in  the  gravitational  geometry. 

Since  we  want  small  changes  in  the  background,  we  must  ask  for  small 
amplitude.  Further,  if  its  frequency  is  low,  its  wavelength  is  large, 
and  in  the  extent  of  one  wavelength,  spacetime  will  be  noticably  curved. 

We  can  hardly  expect  the  signal  to  travel  geodesically  if  it  does  not 
see  spacetime  as  locally  flat,  so  we  probably  should  consider  only  high 
frequencies ^ 67 \  These  are  physical  arguments  designed  to  explain  the 


, 
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fact  that  the  field  equations  in  this  type  of  situation  only  become 

simple  enough  to  solve  easily  after  the  application  of  both  the  weak  field 

and  geometric  optics  approximation. 

Combining  weak  fields  and  geometric  optics  unfortunately  appears  to 

reduce  mathematical  rigour.  We  will  expand  all  quantities  into  double 

series  in  independent  parameters  Aand  6.  Much  less  is  known  about 

the  convergence  of  such  double  series  than  about  single  parameter  series . 

Furthermore,  the  order  of  taking  the  zero  limits  of  the  two  parameters 

can  be  nontrivial.  Classically,  the  energy  density  of  an  electromagnetic 

wave  depends  on  the  amplitude  only,  so  to  drive  energy  density  to  zero  we 

must  take  A  to  zero,  regardless  of  £.  Further,  the  expressions  found 

here  tend  to  have  unexpected  properties  when  we  take  £  to  zero  while  A 

//o) 

is  finite,  so  we  will  assume  that  A  is  not  much  greater  than  » 


A.  <  i 
£  ~ 


symbolically: 

(5.10) 

Consider  now  small,  high  frequency  electromagnetic  and  gravitational 
perturbations  of  a  purely  gravitational  geometry.  We  already  have  applied 
the  weak  field  approximation  and  found  eq. (4.19)-(4.22) .  For  high 
frequency  within  this  we  take: 

(5.H) 
(5.12) 

h  is  slowly  varying,  of  order  unity,  and  we  have  not  written  out 
ab 

1 

possible  terms  in  g  which  are  higher  order  in  A.  Similarly. 

F  =  f  exp(^y)  +  cc  +  . . .  .  13) 

ab  ab  r  c 

Recall  that  the  two  types  of  perturbation  are  independent  here:  we  can 
consider  them  separately,  and  even  go  back  and  use  different  s  and  £for 


gab  =  gab  +  Kb  +  0( 

®ab  =  (hab  exp(f }  +  °C)  +  0{£)' 


each. 


. 
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The  first  order  term  in  the  perturbed  connection  XP  a  is  now: 

be 


^ - 1 0  »*  <t>  ..  *  v 

^bc  kd^  +  CCJ  +  *  *  * 


(5.14) 


x 

e 


fc  exp  +  cc 


+ . . . 


The  leading  term  is  of  order  which  we  are  tempted  to  declare  as  unity 

\ 

Let  us  keep  this  —  term  separate  from  the  order  unity  terms,  even  if 
we  expect  they  may  be  of  similar  size.  Then: 

exp  (f)  i  or°b  k0  -y^0  V  +  oo]  +  --- 

exp  (^)  (-1)  (I  °d  hdb  ko  ka  +  §  c  ’ 

+  ccJ  + . . . 


R  -X 

ab  e2 


82L 


-  h 


ab 


o  cdT  T  o  cd  , 
g  kc  kd  -  g  h 


h  k  k, 
ac  d  o 


cd  ka  hi-* 


(5-15) 


There  is  no  —  term  in  T^,  so  the  term  displayed  here  must  be  zero 


This  is  easily  solved  by  k  null  with  respect  to  the  background: 
Icdkokd  =  0, 


(5.16) 


and  the  "polarization  tensor"  ,  h  traceless  and  orthogonal  to  the  prop¬ 
agation  vector  (with  respect  to  g^): 


o  cd  ,  n 
g  h . ,  =  0 


cd 


o  cd 
g  h 


k  =  0 
ac  d 


(5.17) 

(5.18) 


Then,  by  the  same  argument  used  before,  k  is  also  geodesic.  These 

cl 

are  the  standard  results  of  the  geometric  optics  of  gravitational  waves. 
Considering  the  electromagnetic  case  separately,  the  curl  equation 


is  solved  by: 


=  exp  (^p)  (  *°a)  +  cc  +  ... 


Fab  =  exP  (¥)  w[a  V)  +  00  ’ 

and  the  divergence  equation  is  then: 

[a  'S]  V*  “  +'-- 


/ 1 s \  /i  o  3/C  o  bd  /  . 
0  =  exp  hr)  (?  g  g  “>■ 


(5-19) 

(5.20) 

(5-21) 
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implying: 

o  ab 


*b 


°  ab  y  n 
g  «  k^  =  0 


=  0  (5.22) 

(5-23) 

The  propagation  vector  is  null,  geodesic,  and  orthogonal  to  the  polariza¬ 
tion  vector,  with  respect  to  the  background. 

At  this  point  we  should  comment  further  on  the  relative  size  of  "X  and 
£ .  The  energy  density  of  the  electromagnetic  wave  is  of  order  (the 
leading  terms  in  T  ,  ie, small,  but  the  energy  density  of  the  gravitational 
wave  is  difficult  to  specify  .  Note  that  the  contracted  curvature  and  the 

stress-energy-momentum  are  identified  in  Einstein's  equations,  yet  a 

a  \ 

typical  component  of  R  is  of  order  which  may  diverge  unless 

£  2 

we  make  a  stronger  assumption  than  eq.(5*10),  or  rescale  h^— »£  h^, 


ab1 


so  that  R  ,  ^X(ln  this  context ,  the  rescaling  is  a  mathematical  trick, 

cl  D 

and  does  not  alter  any  results  found  so  far) .  Rigorously,  the  relative 

X 

acceleration  of  nearby  test  particles  will  diverge  if  —  is  of  order 
unity  and  R  It  appears  we  must  be  careful  with  the  relative 

ci  D  gc 

size  of  our  two  parameters  to  get  well-behaved  gravitational  radiation. 
-General  Background 

When  there  is  a  background  F  ^  the  mathematical  complication  increases 

V  f 

We  apply  high  frequency  assumptions  for  g^,  F  ^  to  the  general  back- 

i 

ground  perturbation  equations,  eq . (4.28)-(4.33)  (starting  with  to 
automatically  solve  the  curl  equation) : 

(5.24) 

*  £  a' 

'is 


A  =  exp  (f)  Cf  «  )  +  oe  +... 


Fab  =  6XP  (fW[a  V]  +  CC 


+  ... 


(5.25) 

(5.26) 

We  have  chosen  the  particular  forms  here  on  energy  density  grounds.  To 
mate  the  equations  look  reasonable ,  we  ask  that  the  leading  terms  in  energy 


gab  =  exp  (^r)  (-  £  hab)  +  °=  +  •  •  •  • 


. 


■ 
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density (using  the  term  loosely) ,  which  are  simply  the  leading  terms  in 


Einstein's  equation,  be  of  order  \  if  to  and  h  ^  are  of  order  unity. 

» 

This  is  mainly  a  mathematical  trick  to  avoid  the  problem  of  divergent  R  ^ . 
We  can  let"X~ £n,  and  find  equivalent  equations  for  an  expansion  in  "X', 
with  all  quantities  rescaled  by  appropriate  powers  of  £  .  In  the  strong 
field  high  frequency  case  discussed  earlier  such  rescaling  is  not  possible . 

I  » 

Much  more  important  here  is  the  relative  size  of  g^  and  F^ .  Physically , 

we  can  expect  three  distinguishable  cases.  If  the  source  is  primarily 

electromagnetic  we  may  expect  the  resultant  perturbation  to  have  stronger 

(in  energy)  electromagnetic  than  gravitational  part,  and  conversely,  a 

neutral  massive  source  might  generate  a  mainly-  gravitational  perturbation. 

Finally,  a  charged  massive  source  likely  generates  the  general  case.  ■  - 

These  arguments  are  little  more  than  extrapolations  from  simpler  cases, 

but  if  we  find  freedom  to  adjust  relative  strengths,  it  will  be  tempting 

to  say  the  physical  cases  are  appearing  .  Since  geometric  optics  always 

decouples  the  fields  from  the  sources,  the  question  cannot  be  settled  in 

our  framework.  As  far  as  working  with  eq. (5.24)-(5*26) ,  let  us 

consider  that  and  h  ,  may  be  as  large  as  order  unity,  but  they  may  be 

a  ab 

much  smaller  (so  that  the  ratio  of  their  magnitudes  may  vary). 

Let  us  begin  by  assuming  c o  is  order  unity.  We  find,  from  eq.(4.28): 

cl 

°ef 


Tab  =  exp(^)  +  F\<‘kko3'+  *  ®ab  ^  “teVp 

+  cc  +  ... 

(again,  indices  of  Fab  have  been  raised  with  ga; .  Meanwhile: 


(5-2?) 


ffcG  =  I  §ad  exp(f)  (-16)  (  hbdkc  +  bdA  "  W  +  cc 


Rab  =  4  g°d  exp(^)  (hadVc  +  hdbkakc  "  habkdkc 

-  hn  k  k,  )  +  cc  +  ...  . 
dc  a  d 

For  the  electromagnetic  divergence  equation  (eq(4.32)): 


(5.28) 

(5.29) 


. 

■ 
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Fab  = 


oac  obd  ' 
g  g  F 


cd 


'ab  _  /  oac  obd  '  n 
F  ;b  -  <  S  S  Fcd  )  + 


(5.30) 

(5.31) 


oac  obd  /is\ /i\ /  _  T 

=  g  g  exp(-g-)  (?)  ( 

=  0  . 


"clVS  )  +  00  +  •  •  • 


This  is  the  same  form  seen  in  the  case  of  no  electromagnetic 
background: 

§M  kb  kd  =  0 

ibd^-d  =  °  • 

The  leading  terms  of  Einstein's  equation  then  imply: 

•|gC^(  h  ,k,  k  +  h,,  k  k  -  h  k  1  )  =  F  cu>  k  +  Fc  u>  k 

ad  o  c  db  a  c  dc  a  d  '  a  [jc  oj  b  [a  cj 

0  fief  T 

+  sabF  wEekf]  ' 

I 

It  is  possible  that  the  leading  terms  of  T  ^  (eq(5.27))  are  zero, 

'  o 

but  this  requires  special  coupling  of  F  ^  to  Fq>i,  eg: 


(5.32) 

(5.33) 


(5.34) 


ab! 


obc  2t 
g  F  ,  k  =  0  . 

ab  c 


(5.35) 


This  seems  restrictive.  For  instance,  no  radially  travelling  high 


frequency  F  ^  in  the  Reissner-Nordstrom  geometry  can  satisfy  eq(5*35)* 

Since  it  appears  that  nonzero  eq(5.27)  is  the  general  case,  we  find 

that  h  ,  must  be  order  unity,  and  that  it  must  satisfy  a  rather 
ab 

complicated  equation  linking  it  to  the  electromagnetic  quantities. 
Since  the  electromagnetic  perturbation  travels  as  freely  as  in  the 
no-electromagnetic-background  case,  but  the  gravitational  perturbation 

V 

cannot  be  set  to  zero  because  it  is  coupled  to  F  an  appealing 
interpretation  is  that  this  configuration  is  primarily  an  electro¬ 
magnetic  disturbance,  generating  a  sympathetic  gravitational  distur¬ 
bance  through  interaction  with  the  background. 


Now  let  us  find  the  minimal  F ^  required  when  there  is  a  gravita- 
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tional  perturbation,  due  to  interaction  with  F  We  take  eq(5.26) 
with  h^  of  order  unity,  so  that  R  ^  is  of  order  ,  Obviously,  if 
co  is  any  smaller  than  order  unity,  the  leading  terms  in  Einstein's 

cl 

equation  are  purely  gravitational,  signalling  a  gravitational  wave 
not  interacting  with  any  electromagnetic  field,  and  we  find: 


ocd 

h  ,  =  ° 
cd 

g 

ocd 

g 

h  k  =  0 

ac  d 

ocd 

tv 

0 

tv 

11 

0 

■ 

g 

(5.36) 

(5.37) 

(5.38) 


'ab 

In  the  electromagnetic  equations,  F  contains  terms  like: 


oac  ,bd  °  oac  obe  odf  , 

g  t  Fcd  =  -g  g  g  k 


ef  Fcd  ’ 


'ah  ?  '  '  ab 

so  F  is  of  order  fc  even  if  there  is  no  F  ,  .  Then  F  , 

ab  ;b 

I 

terms  in  order  £  .  If  F  ,  is  zero: 

ab 

'ab  oac  ,bd  ft  ,  ,ac  obd  ft 

F  =  g  k  Fcd  +  k  S  Fcd  +  . . . 


(5-39) 
has  several 


ftb/  oae  ocf 


=  §df  sef )  +  -§ 


4ef  )  + 


(5.^) 


/isNC2/  gaf  obe 
=  exp  (-g- )  £  (  F  g  h 

and  eq(4.32)  becomes: 


ofb  oae  ,  \  , 

ef  +  F  s  hef  )  +  00  +  . . .  , 


n  _  p&b  ,  §ac  p  t>  +  n  a  + 

0  "  F  ,b  F  I  cb  +  *  I  cb  '*•  ’ 


(5.41) 


which  simplifies  using  eq(5*36)-(5.38)  to: 

0  =  |ae  hef  Fbf  k^ 

(leading  term) .  This  is  an  extra  restriction  on  the  gravitational 

I 

perturbation  (which  is  undesirable)  that  can  be  avoided  if  F^  is 
2 

of  order  £  : 


(5-42) 


'  /isx  _2 

Fab  =  exp(r )  £ 

'ab 


>  k,.,  +  cc  +  ...  (5.43) 

aD  “  ■  c.  -  La  oj 

J*>  =  exp(f  )e2(  r°  +  h=a<  Fad  +  §ac  Fdb)) 


+  CC  + 


0  •  • 


(5.44) 
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Instead  of  eq(5.^2)  we  find: 

_  i  oac  obd  T  _  oac  °bd  /  s 

0=TS  g  %  kc  -  g  hod  F  ,  (5.45) 

I 

indicating  that  F  ^  is  coupled  to  the  background.  From  the  same  kind 
of  argument  as  before,  this  appears  to  be  a  source-generated  gravita¬ 
tional  perturbation,  producing  an  electromagnetic  disturbance  through 
coupling  to  the  background. 

I 

Between  the  two  cases  discussed  so  far  lies  F  ,  of  order  E,  h  , 

ab  ab 

of  order  unity.  On  the  basis  of  our  physical  arguments,  this  should 
mean  both  types  of  radiation  are  driven  by  active  sources,  in  which 


case  the  £  appearing  in  the  energy  density  of  F  ^  may  seem  strange . 

It  simply  means  that  the  relative  amplitudes  of  the  two  disturbances 
for  which  geometric  optics  can  identify  both  as  source-driven  (if  this 
is  the  proper  interpretation  of  this  case)  are  frequency  dependent. 
Inspection  of  the  field  equations  for  this  case  reveals  that,  in  the 
leading  order,  the  two  radiation  fields  are  independent  of  each  other, 
and  propagate  without  explicit  regard  for  the  background  electromagnetic 
field  (ie,  along  null  geodesics  as  found  in  the  case  of  a  purely 
gravitational  background) . 

In  the  case  of  independent  leading  orders  (the  last  case  discussed) , 
we  can  again  use  independent  s^ ,  s^,  and  ^ for  the  fields.  In  the 
other  two  cases,  we  have  assumed  that  the  sympathetic  field  (generated 
by  the  interaction  of  the  applied  field  with  the  background)  is  of  the 
same  frequency  as  the  applied  field,  and  have  written  the  equations  for 
the  sympathetic  field  with  that  assumption.  We  have  not  solved  those 
equations  because  they  involve  the  background  and  must  be  considered  on 
a  case-by-case  basis.  If  in  some  particular  case  the  leading  order 
equations  in  the  sympathetic  field  do  not  have  reasonable  solution,  then 


- 


the  frequency  assumption,  and  the  use  of  geometric  optics,  should  be 
reconsidered  for  that  particular  case. 
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In  considering  all  three  cases  we  see  that  a  high  frequency 
perturbation  which  is  applied  to  the  background,  and  not  generated  by 
the  interaction  of  another  field  with  the  background,  travels  along  null 
geodesics  of  the  background,  in  first  approximation. 

Geometric  Optics  In  BMB 
-  Minkowski  Background 

We  noted  that  the  EMT  equations  for  the  first  order  of  a  perturbation 
of  Minkowski  space  have  plane  wave  solutions,  without  any  restriction  on 
frequency.  Mimicing  a  standard  method  of  solution  in  EMT  for  the  analo¬ 
gous  case  in  BMB,  we  ask  for  harmonic  coordinates: 


g 


ab  nc 


r 


ab 


=  0 


(5.46) 


at  least  in  first  order,  where  this  implies: 


ab  »  i  „ab  » 

V.  ^(bc),a  2  U  ^(ab),c 


(5.47) 

Then,  note  that  if  the  first  order  BMB  fields  satisfy  eq.(5.4?)  and 
(using  eq(2.11) , (2.17)  for  the  generalized  electromagnetic  field  and 
potential) : 


f  1  • 

E  _  =  A  .  -  A, 

(5.48) 

ab  a ,  b  b ,  a 

2  1 

rr  a  =  o 

(5-49) 

a 

ab  I  n 

?  Aa,b  =  c 

(5-50) 

2  •  A 

O  *(ab)  =  0 

(5.51) 

(C  is  a  constant)  then  they  satisfy  eq(4.5l)_(4.53) •  Eq (5*47)  (5*51)  are 
(setting  G  to  zero)  just  the  usual  EMT  equations  in  Lorentz  gauge  for 
these  circumstances,  and  they  have  the  standard  plane  wave  solutions: 


A  (x)  =  co  exp(ik  xc)  +  cc 
a  a  c 


(5.52) 


. 

'  '  ■  i  i  i  Hli 


■ 


84 


gab(x) 


=  hab  exp(ipcx°)  +  co 


(5- 53) 


where : 


T?"13  ka  ^  =  0 

?ab  pa  pb  =  0 

rf0  k  w  =  o 

t  a  b 

ab  u  _  i  ^ab  , 

?  pa  bbc  -  ^  U  hat  p< 


(5-54) 

(5.55) 

(5-56) 

(5-5?) 


and  k  ,  p  ,  a>,  ,  h,  are  constants,  k  and  p  define  straight  lines’ (of 

cL  ct  U  UL  3,  3/ 

propagation) ,  which  are  the  geodesics  of  Minkowski  space  (to  first  order 

k  ,  is  itself  zero) . 
a;b  7 

In  chapter  4  we  made  some  comments  about  the  second  order  of 

perturbation.  In  particular,  if  the  first  order  is  purely  electromagnetic 

(ie,  g^^)  is  zero),  then  the  seconds  order  can  be  purely  gravitational 

(ie,  g_  zero),  satisfying  EMT's  Einstein  equation  with  correction 
l_aDj 

terms  of  order  k  (This  k  is  the  new  coupling  constant) .  In  the  case  of 

first  order  electromagnetic  plane  wave,  eq(5*52) , (5*5^) » (5*56) ,  the 
2 

order  k  terms  in  the  second  order  gravitational  equation  are  all  zero 
(see  Appendix),  so  that  only  Einstein's  equation  is  lefts 

(5*58) 


I V  I! 

R  ,  /  ,  \  =  8ttT  ,  /  ,  \  . 

ab(emt)  ab^emt; 


Thus,  for  plane  electromagnetic  waves  perturbing  Minkowski  background, 
EMT  and  BMB  are  not  distinguishable  to  second  order  in  amplitude. 

-  Gravitational  Background 

We  noted  in  chapter  4  that  in  a  purely  gravitational  background  the 
BMB  g^  ^  perturbation  is  independent  of  g »  anc^  sa^as-^aes  equations 
identical  to  those  of  the  analogous  general  relativity  case.  The  results 
for  radiation  in  this  case  are  well  known v  ' ,  and  have  already  been 
discussed.  For  the  electromagnetic  field, 


* 


' 


85 


S —  exp(“e)  (“S  ^[abQ^  c<^  "**  "  ’  ’  (5*59) 

we  must  reconsider  our  arguments  about  relative  amplitudes.  In  EMT  we 

*  2 

argued  that  in  high  frequency  g^^N  must  be  no  larger  than  0(  £  )  and 

f 

F  ^  no  larger  than  0(1)  so  that  energy  as  it  appears  in  Einstein's  equation 
is  0OO.  Here,  g^^  and  are  on  an  entirely  equal  footing,  and  if 

we  still  ask  that  typical  curvature  components  have  amplitude  ,  then 

'  2 

g^^-j  must  be  of  order  £  .  This  appears  to  be  a  decided  change  from  EMT, 
but  it  remains  to  be  seen  whether  it  generates  new  properties. 

Let  us  now  examine  the  leading  order  terms  in  the  field  equations: 

1 

a 


r  be  =  -P  <¥>  <-  §  ad(hCbd]kc  +  V]1^ 

+  h[>c]kd')  +  cc  +  •  •  •  » 
but  this  must  be  traceless,  implying: 
o  cd 


g 


kc  h[M]  =  0 


(5.60) 


(5.61) 


Now: 


R[at]  '  r:ab],o  + 

=  exp  (^r)  j  g  d(h|2adJkb  +  h [db] ka  +  h[ab]k4  ke 


(5.62) 


+  CC  +  ... 

/is\  1  o  cd  T  T  ■,  ,  , 

=  exp  (-g-)  2  g  k0kdh[ab]  +  00  ''■••• 

while  I  is  of  order  £2.  Therefore  the  field  equation  on  the  curvature 
ab 


is: 


O  Cd  -  T  1  r\ 

g  kckdhCab,e]  ~ 

which  is  solved  by  : 

o  cd  T  T  n 
g  k  k  =  0 
c  d 


(5.63) 


(5.64) 


or 


hr  ■>  =  0 

[|ab ,  e  J 


(5.65) 


If  eqi.  1(15.  65)  holds  theii  there,  exists  some  vector  v a>  not  proport 


.. 
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ional  to  k  ,  such  that: 

cl 


h[ah]  V[akb]  ’ 


(5.66) 


"but  we  must  have: 


(5.6?) 


which  can  only  he  solved  l?y; 


(5.68) 


and: 


(5-69) 


The  conclusion  that  can  now  he  drawn  is  that  to  satisfy  eq.(5.63)  k 


a 


must  he  null  with  respect  to  the  background.  Then  hy  working  through  the 
same  argument  used  in  EMT,  we  find: 


(5.70) 


that  is,  the  paths  generated  hy  k  are  geodesics  of  the  background,  to 
first  order. 

Notice  in  the  above  argument  that  we  have  not  yet  seen  anything 
to  do  with  the  noncompatihle  BMB  connection.  Since  the  background  is 
purely  gravitational,  it  is  ruled  by  general  relativity,  where  the 
connection  is  compatible  with  the  metric.  To  first  order,  only  the  back  - 
ground  connection  appears  i>n  the  geodesic  equation.  To  gain  information 
on  whether  radiation  does  follow  the  BMB  noncompatible  connection  we 
must  either  solve  higher  orders  here  (which  would  require  detailed 
knowledge  of  the  background)  or  leading  order  in  the  general-background 
case,  which  we  will  attempt  to  do. 

-  General  Background 

Using  the  high  frequency  form: 


(5.71) 


. 


8? 


(with  both  symmetric  and  skew  parts)  in  eq.  (4. 65)-(4. 70 ) ,  we  find  for 


the  connection  (due  to  the  high  frequency) : 

o  1  d  '  1  ' 

g(cd)  ^  (ab)  2  ^g(cb),a  +  g(ac),b  g(ab),c^ 


o  p  d  o  p 

S[ldbJ  Of).  "  g[adJ  1  [be] 

0  p  d  _  i  /  '  '  \ 

S(cd)  1  [lab]  2  ^Stcb],a  +  g[ac],b  +  SCab],cJ 


n  d 


(5.72) 


__  o  r—»  d  p  d 

g[db]  '  (ca)  gpadl|  (be) 


'n  d 


(5.73) 


rSL  o  s»"b 

we  must  contract  these  with  if  ,  the  inverse 

of  g^-^ ,  which  does  exist  for  a  well  behaved  background  metric.  However, 

doing  so  resembles  "raising  indices  with  the  symmetric  part  of  the  metric", 

contrary  to  the  spirit  of  chapter  3*  Unfortunately,  we  have  no  choice  if 


we  wish  to  go  on.  It  does  at  least  seem  consistent  to  assume: 

r^c =  exp(|?)  (- 18  +  cc  +  •  •  •  • 

Then: 


(5.74) 


w  c  o  ce 

^ (ab)  ^ 


^h(eb)ka  +  h(ae)kb  '  h(ab)ke} 


o  yd  o  yd 

‘  s[db]  °l&f\  s[ad]  "[.bej. 


(5.75) 


c  o  ce 
*  [abl  =  r 


[^hCeb3ka  +  h[ae]kb  +  hr-'--tkJ 


'[ab^[  e' 


o  v  d.  o  Vd 

"  SL^1  d  (ea)  g[ad]  d(be)_ 


(5-76) 


These  forms  do  not  allow  us  to  immediately  express  in  terms  of  h  ^ 
and  k^  because  of  the  interaction  of  with  through  • 
Proceeding  on  to  the  other  equations: 


°  ”  6>c] 

=  ^°eh[ae3ko 


o  v  d  o  ce  o 
g[dc3  (ea)  ^ 


Y  d 

[ad]  °(ce) 


(5-77) 


’  involving  no  dsxivativos ,  is  much  small@x  "than  "th©  loading  cunvatuxs 
ab’ 

terms,  and  can  be  ignored.  Thus  the  curvature  equations  (leading  terms) 


are : 


•  '  ■  -  ^ 
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0  =  ^*0  -  *<  A  +  XfaK) 

^[a.b f  eXp(T}  ko  ^Eab]  +  00  +  ••• 
With  the  reasonable  assumption: 

'  is 

=  exp(-g-)  (-ieab)  +  cc  +  ... 


we  find: 

2  a  p,  k  —  k ,  p ,  -t  | 
[b  c]  d"[bo] 

which  is  solved  by  : 

%  =  kd  pd  qb 
*[bc]  =  2  pd  q[bko]  • 


(5-78) 

(5-79) 

(5-80) 

(5.81) 

(5.82) 

(5.83) 


Then  the  trace  equation  says  k  and  p  are  orthogonal,  which  means: 

cl 


k  yrc ..  =  o 

c  [abJ 


(5.84) 

[abj'S  -*-eac^-n&  order  is  itself  zero.  Eq.(5.84)  is  probably  the 


and  Rj 

strongest  thing  that  can  be  concluded  (rather  than  any  of  the  intermediate 
steps).  Using  eq.(5*76): 


c  ~  Qab] 
i  /,  o  ce 


0  =  k, 

(hCeb]  « 

,  .  o  ce  T  T  \  o  v/  a  o  ce  T 

+  h[ab]  V  koke  -  SCdb]  ^(ea)  «  k< 


k  k  +  hr 
c  a  LaeJ 


oce,  T 
*  kokb 


(5-85) 


o  n/  d  o  ce. 

®[ad^J  (be)  c 

At  this  point,  further  progress  becomes  difficult.  There  is  no  obvious 
solution  to  eq.(5.85),  but  we  are  hoping  for  simple  results,  perhaps: 

0  :  (5-86) 

(5.87) 

[aej  c 

and 


oce.  , 
n-'  k  k 
^  c  e 


oce,  T  _  n 
^  h[ae]ko  -  0 


o  ce 


Y/  Nk  =  0  . 
u  (ea)  c 


(5.88) 


These  three  together  will  solve  eq.(5*85)»  e(l*(5*86)  in  eq.(5*7?) 

implies  a  coupling  of  first  order  to  the  background  and  eq.(5.88)  and 


/ 
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(•5-75)  together  imply  a  different  equation  coupling  first  order  to  back¬ 
ground.  Such  (extra)  couplings  make  the  simple  results  we  are  looking  for 
unlikely . 

The  same  simple-minded  techniques  can  be  applied  to  the  equation  in 
symmetric  curvature.  To  satisfy  eq.(5.78)  we  might  assume: 

Y  (ac)  =  0 


kc  ^(ab)  =  0 
from  which  we  find: 

o  =  y  ® 

(ac) 

1  oce,  T  oceo  y  d 

2  r  h(ec)ka  ^  s[dc]  ®  [ea] 


0  =  kc  *(ab) 

i  /,  o  ce.  T  ,  ,  O 

=  2  (h(eb)  ^  koka  +  h(ae)  * 


(5.89) 

(5-90) 


(5-91) 


o  ce 


kckb 


(5.92) 


,  oce.  ,  \  o  yd  o  ce, 

(ab)  ^  c  e  ^[cLb]  [ea]  ^  c 


(5-92) 

In  the  same  way  as  before,  trying  to  assume  simple  forms  here  results 


o  V  d  o  ce, 

'  S[ad]  IbeH't'  kc 


in  extra  coupling  of  background  to  foreground.  We  fall  into  calculational 
confusion. 

Having  been  frustrated  in  this  approach,  we  can  try  changing  the 

V  f 

relative  size  of  and  g£ab3  (in  analo&y  with  the  EMT  section). 

If  we  take: 


[abj 

exp(^)  (-£2h|-ab:1)  +  cc  +  ... 

(5-93) 

(ab) 

0  (  S3) 

(5-94) 

V  c 

o  ce  o  yd  oce  o  yd 

-  "  r  s[db]  "  tea]  r  SCad]  X[be: 

(5-95) 

3  (ab) 

y  c 

0  Lab] 

=  r°d  j(h|-db-]  \  +  hjad]  + 

(5-96) 

then: 


!  (.tf~  iY  <m\a)  f 
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ref(  ltfb3lrjd]^a:  +  S[fb]Sl;aj]y[:deD  +  ®[af]S[jd]yCeb] 


o  o 
+  S  r^-PiS 


c 

Here,  instead  of  expressed  in  terms  of  hj-^-j  and  ka  (which  we 

expected)  ,  we  find  a  complicated  expression  in  ^j“a^  and  'the  background 

electromagnetic  field.  Without  an  easily  solved  relation  between  g  and 

ab 

^ab>  it  becomes  difficult  to  draw  any  conclusions  from  the  rest  of  the 
field  equations .  The  author  has  had  little  success  in  assuming  simple 
forms  (in  analogy  with  EMT)  in  this  case  or  in  the  case  of  g(ab)  larger 


than  g 


Gaby 


In  summarizing,  the  failure  of  this  approach  we  can  make  several  remarks . 

Whereas  in  the  earlier  cases ,  there  seemed  to  be  a  fair  analogy  between 

EMT  and  BMB  radiation,  there  is  no  similarity  apparent  here.  The  elementary 

techniques  we  have  used  seemed  to  work  with  ridiculous  ease  in  EMT  by 

contrast  with  this  section.  We  found  that  when  there  is  a  leading  order 

g  _  it  does  not  generate  any  generalized  electromagnetic  potential  at 
Lao  J 

f 

the  order  expected,  so  the  leading  order  S|-ab-j  disappears  in  the  EMT  limit, 
probably  placing  this  disturbance  in  the  "special"  category  (see  eq(2.96), 
(2.97))  discussed  in  chapter  2.  Since  we  have  as  much  trouble  with  S(ab) 

as  with  gr  , -..this  is  little  consolation.  Lastly,  note  the  persistent 
[abj 

„  oab  ,,  ,  „  oab  o(ab)  T, 

appearance  of  r  »  versus  the  absence  of  g  or  even  g  .It  appears 

from  the  work  done  here  that  /b,a^>  is  more  important  for  solving  the 

equations  than  either  of  its  competitors.  Whether  it  is  more  important 


physically  is  yet  to  be  seen. 


Note  further  that  if  we  had  found  k  null  with  repect  to 

ci 


o  ab 


then  k  would  be  geodesic  with  respect  to  p  / 

3/  ^ 

o  ab 


c 

ab)lc: 


r 


k  k  =  0 
a  c;b 

ole 


(5-97) 


■ 
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and  would  follow  a  path  which  is  an  extremal  of  arc  length  (the  proof 
is  exactly  that  used  in  EMT) .  On  the  other  hand,  if  k  is  null  with 

cl 

respect  to  g  : 


0-Sa\kb 

■  »<‘b,w 


(5.98) 


it  is  null  with  respect  to  g 


o(ab) 


Corresponding  to  g 


(ab) 


there  is  a 


symmetric  connection  (call  it  s  )  defined  by 

.(ab) 


g 


;c 

lc2 


=  0 


(5.99) 


and  this  time  k  follows  geodesics  of  R-f  n,  „  : 

a  10  (bc)lc2 

g(at)kk  .  =  0  . 

a  c;  b 

olc2 


(5.100) 


By  the  same  type  of  calculation  done  to  show  that  Pf  and  Pf  n  ,  are 

1  be  (bc)lc 

not  related  by  a  projective  transformation^^ \  we  can  see  that  r(bc)l=2 

rSL  r-i  cL 

be  anc^  (bc)lc  ’  ^us 

the  geometric  optics  formulation  we  have  developed  here  discriminates 

against  the  most  expected  case:  if  high  frequency  radiation  travels 

along  null  geodesics  of  "  in  the  first  approximation",  then  the 

equations  here  will  be  of  complicated  form.  Possibly  the  standard 

analysis  could  be  modified  to  single  out  null  geodesics  of  instead 

of  Pf  ,  which  the  work  here  seems  to  be  doing  (but  the  form  of 
(bc)lc 

such  modification  is  not  obvious  to  the  author) . 

-  Other  Cases 

The  case  of  high  frequency  without  weak  fields  was  not  presented 
earlier  because  it  contains  difficulties.  The  obvious  form  to  assume  is: 

SCab]  =  exp(¥)  h[ab:  +  00  (5-101) 

as  the  generalized  electromagnetic  field  in  some  region  where  the  BMB 

equations  are  satisfied  exactly .  The  field  equation  relating  g^  and 


« 


. 
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cL 

P^c  then  implies: 


r(ab)  =  ^  U  (s(eb),a  +  «(ae),b  '  «(ab)fe} 
-  (exp  (~)  (hrj_rri.  +  hP..an  Frt  i)  +  cc) 


„ce 


S)  T  nCad3  '[be]' 


(where  '7^'°  is  the  inverse  of  g,  N ) 

r-  tab )' 


(5.102) 


n 


Cab  J 


=  r 


ce 


1  i 

2  £ 


*  gXp(v)  (hCebn  ka  +  VeA 


(5.103) 
+  0  ( 1 )  . 


+  h[ab]ke)  ■  exP(/)  (h[db]  ^(ea)  +  h[ad]  2(be))  +  cc. 

If  we  substitute  eq. (5.102)  into  (5-103),  or  vice  versa  ,  we  find  an 

expression  with  (or  related  to  itself  in  a  complicated  man¬ 

ner.  Now,  we  would  like  to  assume  as  little  as  possible  about  S^y 

cl 

Since  the  curvature  R  ,  ,  involves  second  derivatives  of  e  ,  .  it 

bed  toab 

appears  that  curvature  is  of  order  —  ,  due  to  gj-ab-j  (  we  will  discuss 
the  reasonableness  of  this  later) .  If  we  ask  that  R3^^  be  no  larger 


than  this,  then  g 


(ab) , c 


should  be  no  larger  than  0(-=r) .  If  it  is  of  this 


maximum  size,  then  g 


(ab) ,c 


contributes  directly  to  l"7 


Cbc3  ’ 


Worse 


than  this  are  terms  which  relate  l""^^  to  itself  through  .  In 

these  we  see  factors  like  (exp(-^P))^.  These  indicate  intermodulation, 
where  waves  interfere  to  produce  new  frequencies,  including  higher  and 
lower  frequencies  than  those  of  the  original  waves.  Our  geometric 
optics  forms,  however,  have  assumed  a  narrow  band  of  frequencies,  and 
this  is  not  likely  to  be  the  case  when  there  is  intermodulation.  This 
case  seems  difficult  indeed. 

Possibly  one  of  the  basic  problems  with  this  formulation  of  non-weak 
high  frequency  radiation  is  the  size  of  the  curvature.  Under  normal 
circumstances,  the  curvature  must  be  finite  at  x  for  x  to  be  considered 


part  of  physical  space.  In  the  context  of  high  frequency,  Isaacson 


(70) 


allows  Ra,  ,  to  be  0(-i-)  in  a  particular  situation  in  general  relativity, 
bed  7 

and  then  argues  that  this  will  not  cause  unreasonable  results  in  that 
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particular  situation.  In  the  EMT  non-weak  field  case,  T  ,  is  order 

ab 

unity,  so  the  electromagnetic  field  only  generates  order  unity  curvature. 


The  mathematics  of  that  case  work  nicely  only  if  Ra  ,  is  smaller  than 

bed 

0(^)«  we  assumed  that  fab>c  is  order  unity: 

„  ad  pb  „  db  pa  ,  .  * 

f  ;c  -  ^  ,c  +f  1  dc  +f  rdc'  (5.104) 

and  this  is  easily  so  providing  PaQ  is  order  unity  (ie ,  ~(^)  or 

less)  .  Thus  the  BMB  work  with  R  of  order  -pj  may  not  be  physically 
meaningful.  However,  no  alternate  to  eq.(5.10l)  immediately  presents 
itself  for  the  investigation  of  high  frequency  without  weak  fields. 

So  far  we  have  not  been  able  to  actually  solve  for  paths  of  motion 
when  the  connection  contains  obvious  unified  field  modifications.  As  a 
last  attempt  consider  multiplying  the  field  of  eq.(5.10l)  by  £.  The 
result  will  be  a  one-parameter  weak  high  frequency  field,  which  is  not 
as  well  motivated  as  the  two  parameter  fields  but  is  easier  to  manipulate 


mathematically.  Curvature  will  still  be  divergent,  but  now  only  as 

-j=-  .  Considered  as  a  perturbation  to  a  background  without'  electromagnetic 

field  (and  allowing  a  possible  gravitational  perturbation  of  similar  form) : 

g[ab]  =  eexP^)  h[ab]  +  00  + 

g(ab)  =  ®(ab)  +  eex^¥)  h(ab) 
is  a  solution  in  general  relativity 


+  cc  +  ... 

(71) 


(5.105) 

(5.106) 


'(ab) 


This  case  is  in  a  weak 

(70) 


sense  analogous  to  the  one  considered  by  Isaacson  '  ,  and,  in  a 
different  sense,  can  be  obtained  from  the  earlier  two  parameter 
perturbation  on  a  gravitational  background  by  rescaling  so  that  the 

A-  . 

curvature  is  of  order  2  a-nd  setting  A  equal  to  &  . 

gab  =  o(ab)+  £exp(^I)  (h(ab)  +  h^3^)  +  cc  +  ...  (5-107) 


i — ,  c  i  o(ce)  /O  ,  o  o  n 

r (ab)  =  2  8  (g(eb),a  +  S(ae),b  S(ab),e') 


+  f-exp(^)  §<'Ce'>(h(e'b)ka  +  h(ae)kb  '  h(ab)ke 


k)  +  cc  +  0(£) 


(5.108) 


(5.109) 


I  '  ^  _  f i  o(ce)  r 

I  tab]  -  z  s'  (h[eb]ka  +  h[ae]kb 

+  h[ab]ke')  +  cc  +  0(  S  )  . 

Note  the  symmetric  part  of  the  connection  will  have  order  unity  high 
frequency  terms  if  h^^  is  nonzero,  while  the  torsion  is  order  unity, 
and  high  frequency.  Thus,  P^-^^  ,  which  is  to  order  unity  eq.  (5*108), 
differs  from  the  field  equation  connection  in  order  unity.  Working 
through  the  field  equations  in  the  same  manner  as  was  done  in  the 


two-parameter  gravitational  background  cases 
o  ce 


hr  _,k  =  0 
[acj  e 


CabJ  E  eXp(  £  ^  ^h[abJkekc  +  cc  +  *  *  * 


and  eventually: 
o(ce) 


k  k  =  0  . 
c  e 


(5.110) 

(5.111) 

(5.H2) 


Furthermore,  the  equation  on  is  strictly  in  h^^ ,  so  h^^  and 

^Qab]  are  ^-n(^ePenclen^  1R  leading  order,  we  can  set  h^.^  ^ 0  zero  H  we 
so  desire,  and  we  can  attach  independent  s^  '  and  k'  '  to  the  two 
independent  perturbations.  Anyway,  the  electromagnetic  perturbation 
moves  on  null  geodesics  of  the  connection  compatible  with  the  slowly 


(5*113) 


varying  background  metric: 

S>  c  i  o(ce)  (O  o  o  \ 

r (ab)lc  ■  2®  ^S(eb),a  S(ae),b  S(ab),e^  ' 

At  first  glance,  it  would  seem  that  in  this  case  the  Levi-Civita 
connection  is  definitely  singled  out,  and  the  field  equation  connection, 
being  the  same  size,  is  definitely  avoided.  Unfortunately,  there  is 
another  consideration.  If  we  consider  hp^-j zero ,  h^.^  nonzero, 
we  have  a  gravitational  perturbation  in  general  relativity ,  in  fact ,  the 
leading  order  of  the  case  considered  by  Isaacson.  In  it  we  find  the 
connection,  simply  eq. (5*108),  has  order  unity  high  and  low  frequency 


'  X<u 


. 


, 
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parts,  but  the  radiation  travels  on  null  geodesics  of  the  low  frequency 
part  alone,  eq.(5*H 3) •  Thus  we  once  again  cannot  conclude  anything 
about  the  relative  validity  of  vs  l""'(bC)lC  *  s^-n0e  ^  may  that 

the  low  frequency  part  of  the  connection  (whichever  connection  is  correct) 


is  singled  out  in  this  case. 

Summary 

In  this  chapter  we  have  applied  the  geometric  optics  approximation, 
in  a  form  extended  slightly  beyond  that  of  introductory  relativity 
texts,  to  EMT  and  BMB.  With  it,  radiation  in  EMT  was  found  to  behave  in 
the  intuitively  (even  naively)  expected  manner  in  all  interesting  cases 
tested.  In  BMB,  while  we  did  find  nice  results  in  simple  circumstances, 
we  were  never  able  to  solve  any  case  that  would  have  truly  discriminated 
between  n&oic  and  ,  and  thus  begin  to  resolve  some  of  the 

questions  raised  in  chapter  3*  Possibly  we  could  evaluate  the 
non-weak-field  and  general  background  cases  using  Tonnelat '  s^^'  method 
of  finding  »  but  this  is  beyond  the  scope  of  this  thesis.  We  could 

try  to  test  the  general  background  case  using  the  exact,  spherically 
symmetric,  static  solution  as  a  particular  background.  While  this  might 
help  rule  out  some  of  the  many  possible  radiation  trajectories,  it 
would  be  difficult  to  identify  the  general  rule  from  one  such  example. 
Lastly,  geometric  optics  is  not  the  only  way  to  study  radiation.  The 
study  of  characteristics  (as  mathematical  idealizations  of  shock  wave 
fronts)  has  proven  fruitful  in  general  relativity  ,  but  is  once  again 
beyond  the  scope  of  this  thesis. 


. 
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Chapter  VI 
Concluding  Remarks 
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In  the  past  five  chapters  we  have  progressed  through  a  variety 
of  topics  in  an  effort  to  gain  understanding  of  Bonnor-Moffat-Boal 
unified  field  theory.  Chapter  1  was  a  review  of  the  subject  of  unified 
field  theories  of  gravity  and  electromagnetism,  chapter  2  a  review  of 
BMB  theory  with  a  discussion  of  the  concept  of  duality  and  some  interest¬ 
ing  possible  families  of  solutions.  In  chapter  3  we  tried  to  extend  the 
physical  interpretation  of  all  geometric  quantities,  which  makes 
Einstein-Maxwell  theory  so  aesthetically  pleasing,  to  BMB.  In  chapter  4, 
perturbation  expansions  were  made  for  a  number  of  interesting  cases, 
and  with  them  an  approximate  BMB  magnetic  momopole  was  found.  Lastly, 
in  chapter  5>  we  discussed  BMB  radiation  in  the  high  frequency  (geometric 
optics)  limit. 

Chapter  1  generally  showed  the  troubles  of  unified  field  theories. 
Weyl's  theory  is  an  elegant  gauge  field  extenesion  of  Riemannian  geometry 
but  it  is  difficult  to  find  acceptable  field  equations.  The  five  dimen¬ 
sional  theories  suffer  from  a  lack  of  geometric  interpretation  of  the 
fifth  dimension,  and  they  either  generate  unusual  field  equations  or 
simply  are  equivalent  to  EMT.  Einstein's  nonsymmetric  field  seems  like 
a  logical  generalization  of  general  relativity,  yet  it  has  strange 
exact  solutions  and  equations  of  motion. 

In  chapter  2,  we  saw  that  BMB  theory  is  a  modification  of  Einstein's 
nonsymmetric  theory  with  improved  equations  of  motion.  Furthermore  it  is 
a  generalization  of  EMT,  not  just  of  general  relativity.  It  has  an 
exact  electric  monopole  solution  ,  and  may  not  have  any  magnetic  monopole 
solution  of  similar  form.  To  study  the  possibility  of  a  magnetic  monopole 
we  looked  at  the  concept  of  duality ,  and  found  that  there  is  no  simple 
way  to  formulate  duality  in  a  theory  of  Einsteinian  nonsymmetric  form. 


» 
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Lastly ,  we  pointed  out  several  categories  of  possible  solutions  (eg , 
electromagnetic  field  on  a  flat  background,  or  electromagnetic  field 
with  no  electromagnetic  potential)  that  should  be  investigated  to 
establish  some  of  the  basic  properties  of  BMB  theory. 

In  chapter  3  we  discovered  some  very  basic  problems  in  BMB  theory 
in  that  a  large  portion  of  BMB's  geometric  quantities  and  properties  is 
without  satisfactory  physical  interpretation.  There  are  two  choices  of 
connection,  and  an  inner  product  which  is  not  symmetric  in  its  arguments. 
The  (preferred)  metric  is  not  compatible  with  the  (preferred)  connection, 
so  only  very  limited  families  of  vectors  have  lengths  and  inner  products 
preserved  in  parallel  transport .  This  subsequently  makes  well-behaved 
frames  difficult  to  construct .  The  hermitian  formulation  of  the  theory 
solves  some  of  these  problems,  but  creates  others  of  its  own.  Lastly, 
Killing  vector  theory  is  considerably  complicated  by  the  nonsymmetric 
metric,  and  Killing  vectorfields  appear  not  to  be  associated  with 
constants  of  motion  in  nonsymmetric  field  theories. 

Most  of  chapter  4  consisted  of  perturbation  expansions  of  the  field 
equations  for  various  generic  cases.  The  lone  application  of  these  in  the 
chapter  was  the  derivation  of  an  approximate  magnetic  monopole  solution. 
Unfortunately,  if  there  is  a  corresponding  exact  solution  ,  it  will  likely 
be  of  a  form  difficult  to  find. 

In  chapter  5  we  applied  the  weak  field  approximation  of  chapter  4 
and  the  "geometric  optics"  high  frequency  approximation  to  radiation 
in  five  cases.  For  weak  fields  of  arbitrary  frequency  on  a  flat 
background,  and  weak  high  frequency  fields  on  a  purely  gravitational 
background  we  found  fairly  strong  similarity  between  EMT  and  BMB 
radiation.  However,  for  arbitrary  strength  high  frequency  waves  m 
an  exact  solution,  and  for  weak  high  frequency  waves  in  a  general  back- 
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ground,  we  failed  to  find  simple  linearized  equations,  and  probably  should 
not  attempt  to  draw  any  conclusions.  Finally,  for  a  one  parameter 
weak  high  frequency  case  where  we  allowed  curvature  to  diverge  (within 
limits) ,  we  obtained  simple  results,  but  found  them  difficult  to  interpret. 

The  problems  found  in  chapter  3  limit  the  present  acceptability  of 
BMB  theory .  Until  a  significant  number  of  them  are  solved,  the  theory 
will  not  have  solid  physical  interpretation  of  its  mathematics. 

Without  a  reasonably  complete  physical  interpretation  the  theory  cannot 
be  considered  a  serious  competitor  to  EMT,  and  therefore  it  is  the 
author's  opinion  that  at  present  work  in  the  theory  should  be  directed 
toward  resolving  the  interpretational  problems. 
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M  1^ 

i^k)P  happens  to  "be  (see  Ref. 53)  "the  commutator  of  the  two  vectorfields 


k  and  p  when  they  are  considered  as  differential  operators: 


b  T  a  b 
iP  =  k  p 


(R.2) 
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56.  As  an  example  we  can  take  the  static,  spherically  symmetric  solution, 
solve  for  any  radial  (timelike  or  null)  geodesic,  tangent  k  ,  and 
consider  parallel  propagation  of  a  vector  m  which  is  null  at  a  point 
p  on  the  geodesic .  It  is  not  necessary  to  solve  the  parallel  transport 
equations ,  we  merely  look  at  the  derivative  of  m  ' s  magnitude  at  p  : 


TC/  abN  Tc  ab,Tc  a  bLTc  ab  /D 

k  (sabm  m  >;c  =  k  Sab;cm  m  +  k  gabm  ;Cm  +  k  Sabm  m  ;c.  (R>3) 


Now,  we  know  we  can  parallel  transport  m  along  the  geodesic: 

T  c  a  _  n 
km  =0 
;c 


(R.4) 


but  that  g 


ab;  c 


is  not  zero.  The  term  that  is  left  is  generally  nonzero 


a, 


in  the  case  described,  so  as  soon  as  we  move  away  from  p,  m  s 
magnitude  is  nonzero,  and  it  is  no  longer  null. 
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discusses  the  necessity  of  assuming  the  existence  of  well  behaved 
"steady"  coordinates  in  which  the  fields  have  the  form  we  are  consid¬ 
ering.  For  recent  work  on  more  general  cases  in  EMT  see  Refs.  62,63. 
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64.  There  are  entirely  different  methods  of  studying  radiation.  For 
example,  a  number  of  exact  radiation  field  solutions  are  known  and 
have  been  studied  in  general  relativity  (see  J.  Ehlers  and  W.  Kundt 
in  L.  Witten,  "Gravitation,  an  Introduction  to  Current  Research" 
(Wiley,  New  York,  1962)),  which  are  not  asymptotically  flat,  and 
are  therefore  not  easy  to  relate  to  truly  physical  situations . 

Another  way  to  study  radiation  is  to  study  the  characteristic  or 
singular  surfaces  of  the  theory.  These  submanifolds  represent  shock 
wave  fronts  (see,  for  example,  Ref.  5>  pp.  109  -  112,  284  -  289.  For 
recent  work  in  EMT  see  T»  Tokuoka,  J.  Math.  Phys.  15>  1972  (197^) » 
and  Refs .  therin) .  In  our  methods  we  approximate  physical  radiation 
in  the  sense  that  the  first  order  solutions  we  look  at  approximate  an 
exact  radiative  solution  satisfying  the  conditions  of  the  physical 
situation  (if  such  a  solution  exists) . 

65.  This  argument  is  due  to  M.  Walker,  private  communication. 

66.  See  Ref.  53*  PP*  53  -  59* 

67.  To  be  specific,  one  asks  that  g  be  much  less  than  the  background's 
characteristic  length  scale  of  change,  L,  in  the  region  of  interest. 

We  will  take  L  as  unity,  to  simplify  expressions. 

68.  Some  authors,  eg.  Isaacson,  Ref.  60 ,  assume  a  definite  relation  between 

X  and  £  ,  eliminate  one  of  the  two ,  and  work  with  the  remaining  one 
parameter  expansions. 

69.  Isaacson,  Refs.  60,61,  finds  an  effective  stress-energy-momentum  tensor 
which  is  second  order  in  amplitude  for  first  order  waves. 

70.  See  Refs.  60 ,  6l . 

71.  Or  (more  exotically) ,  in  the  spirit  of  Isaacson  (Refs.  60,6l),  g£a^ 

may  be  generated/modified  by  low  frequency  terms  in  higher  order  cur¬ 
vature  actually  generated  by  the  high  frequency  disturbance.  The 
•  differences  between  these  two  approaches  do  not  appear  in  the  leading 
order  we  are  considering. 


Appendix 


Weak  Fields  in  Minkowski  Background,  Second  Order 

Most  of  the  second  order  quantities  in  the  case  of  BMB  weak  fields 

in  Minkowski  space  were  written  out  in  chapter  4,  including  P  f  (start- 

bc 

ing  1  at  eq.(4.38),  ending  at  eq.(4.5l)).  The  remaining  field  equations  are: 

0  =  <U] 

cd/"  1  '  e  '  '  e  •  i  (A .1) 

=  72  (g[bc],4  '  sCeaJ  T (db)  "  «[bej  T (cd)  -  S(de)  f [bof 


!! 

0  =  R 


(ab)  +  "'"(ab) 


A  d 


=  r(ab)fc  ■  ^ r(ac)lb +  r(L),a)  -  r(:d)  r(‘ob) 

,  n  d  p  C  p  c  pd  ,  8jr  /  cd  ' 

I  (od)  I  (ab)  1  Qad]  Ccbf  fc2  g[da]g[cb] 

'*?ab;?  %  g[cd]g[ef]) 

°  =  C(R[ab]  +  ICab]),o] 

/  n  d  _  8tt  "  \ 

■  [U  [ab],d  k2  S[abf  .c]  1 


(A.  2) 


(A.  3) 


As  it  is  written,  eq.(A.2)  seems  similar  to  Einstein's  equation  in  EMT. 
Recalling  that  in  chapter  4  we  showed  that  first  order  EMT  solutions  on 
Minkowski  background  (Lorentz  gauge)  are  solutions  of  BMB  for  the  same 
case,  using  eq.(2.1l): 


& 


[lab] 


=  kF 


ab 


(F  the  generalized  electromagnetic  field) ,  and  noting  that: 
ab 


p.  a  _  n  a  ^ad  '  p  e  ao.  p  e 

I  (be)  ■  1  (be)emt  '  ”  g  [ec]  1  [db]  V  g[be]  1  [cd] 

we  see  that  eq.(A.2)  can  be  written: 

If  »» 

-  8ttT 


ad  ' 


(A  ,4) 


(ab)emt 


(ab)emt 


-  2Cd(S[;eb]  r [da]  +  g[ae]  H:db]te 


-^°d  (gCe=]Pda]),b-^Cd(« 


r e  ) 

fee]  fdb]  y  ,a 


(A.  5) 
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-  r  c  r  d  -  o 
1  Cadll  1  Cob]  ~  U 

II 

where  is  the  (second  order)  electromagentic  stress-energy-  momentum 

'  ? 

due  to  F  The  extra  terms  are  all  of  order  k  ,  so  this  equation  has  the 

expected  EMT  limit. 

Obviously  the  second  order  equations  involve  coupling  to  first 
order,  so  we  cannot  solve  without  using  particular  first  order  forms.  We 
can,  for  example,  look  at  the  cases  of  only  gravitational  or  only  electro¬ 
magnetic  first  order  perturbation.  These  cases  are  possible  since  the 

* 

perturbations  are  independent  in  first  order.  We  know  that  all  solutions 
of  general  relativity  without  electromagnetic  field  are  also  the  purely 
gravitational  solutions  of  BMB.  From  this,  and  knowledge  of  perturbation 
calculations  in  general  relativity,  we  can  predict  that  a  gravitational 
perturbation  to  a  geometry  without  electromagnetic  background  will  not 
generate  an  electomagnetic  perturbation  in  any  order.  This  is  supported 
by  the  BMB  calculations  on  a  purely  gravitational  background,  and  it  can 
be  seen  here  that  if  g^-j  is  zero,  then  the  equations  in  g^j  have 

V 

trivial  solution,  and  those  in  S(ab)  are  those  of  ordinaiy  general  relativity 

t  n 

for  the  same  circumstances.  In  the  other  case,  ie,  6(ab)  zero,  gj-^-j  can 
again  be  zero,  and  S(ab)  is  governed  by  eq.(A.5). 

Finally,  if  the  first  order  solution  is  an  electromagnetic  plane 

wave,  eq. (5»52) , (5*5^) > (5*56) : 

A  =■§■(«  exp  (ik  x°)  +  cc) 
a  a  c 

?Zab  k  k  =0 
a  b 

^ab  k_  ^  =  0  , 


^a  b 


then: 


=  1  Vb]exp  (1VC)  + 


cc 


(A.  6) 


. 

, 
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Hlata  =  exP(ikexe)  (-?°dkd^Cakb])  +  °0' 

'  A  c 

Ij^abQ  ^as  a  ra^er  complicated,  form,  "but  depends  upon: 


(A.  7) 


(A.  8) 


and  similar  expressions  with  cl>  • s  replacing  some  or  all  the  co  ' s .  These 

^  a  .a 

are  all  easily  seen  to  be  zero,  so': 

gtfc]  ^Cab:  =  °.  (A. 9) 

Similarly : 

^[adl  ^Cob3  =  0  ’  (A,1°) 

and  for  first  order  plane  electromagnetic  perturbation  on  Minkowski 

background  (in  BMB)  ,  the  second  order  gravitational  field  is  governed  by 

EMT's  Einstein  equation  for  the  same  situation. 


' 


) 


